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CONVERSION  FACTORS,  INCH-POUND  TO  METRIC  (SI) 
UNITS  OF  MEASUREMENT 


Inch-pound  units  of  measurement  used  In  this  report  can  be  converted  to 


metric  (SI)  units  as  follows: 

Multiply 

By  _ 

To  Obtain 

feet 

0.3048 

metres 

i nches 

2.54 

centimetres 

kips  (1000  lb  force) 

4.448222 

ki lonewtons 

kips  (force)  per  foot 

14.5939042 

ki lonewtons  per  metre 

pounds  (force) 

4.448222 

newtons 

pound  (force)- inches 

0.  1 129848 

newton-metres 

pounds  (force)  per  square  inch 

6.894757 

ki lopasca Is 

square  inches 

6.4516 

square  centimetres 
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USER* S  GUIDE:  COMPUTER  PROGRAM  FOR 
ANALYSIS  OF  BEAM-COLUMN  STRUCTURES 
WITH  NONLINEAR  SUPPORTS  (CBEAMC) 


PART  I:  INTRODUCTION 


General 


1.  This  report  documents  a  computer  program — CBEAMC — for  analysis  of 
general  beam-column  structures  supported  and/or  loaded  by  components  which 
interact  with  the  displacements  of  the  beam  and/or  column.* 

Organization  of  Report 

2.  The  remainder  of  this  report  is  divided  into  the  following  parts: 

<a.  Part  II:  Describes  the  general  beam-column  system  considered 
and  the  mathematical  model  used  for  analysis. 

b^.  Part  III:  Presents  the  force-displacement  relationships  for  the 
mathematical  model  and  describes  the  computational  procedure  used 
for  solution. 

c.  Part  IV:  Provides  example  solutions  obtained  with  the  program. 

Disclaimer 

3.  The  computer  program  described  in  this  report  has  been  checked  to 
ensure  that  the  results  are  accurate  within  the  limitations  of  the  procedures 
employed.  However,  there  may  be  unusual  situations  which  were  not  anticipated 
which  may  cause  the  program  to  produce  questionable  results.  It  is  the  respon¬ 
sibility  of  the  user  to  judge  the  validity  of  the  results.  No  responsibility 
is  assumed  by  the  author  for  the  design  or  behavior  of  any  structure 

based  on  results  obtained  with  the  program. 


*  CBEAMC  is  designated  X0060  in  the  Conversationally  Oriented  Real-Time  Pro¬ 
gram-Generating  System  (CORPS)  library.  Three  sheets  entitled  "PROGRAM 
INFORMATION"  have  been  hand- i nserted  inside  the  front  cover.  The  present 
general  information  on  CBEAMC  and  describe  how  it  can  be  accessed.  If  pro¬ 
cedures  used  to  access  CORPS  programs  should  change,  recipients  of  this  re¬ 
port  will  be  furnished  a  revised  version  of  the  "PROGRAM  INFORMATION." 
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PART  II:  BEAM-COLUMN  SYSTEM 


General 


i  . 


4.  The  general  beam-column  system  considered  for  analysis  is  shown 
in  Figure  1.  Characteristics  and  assumptions  employed  for  each  of  the 
system  components  are  described  below. 

Beam 


5.  The  following  assumptions  of  conventional  beam  bending  theory 
are  employed: 

£.  The  beam  is  composed  of  a  linearly  elastic,  piecewise  homo¬ 
geneous  material  with  modulus  of  elasticity  E. 

The  beam  is  essentially  piecewise  prismatic,  i.e.,  the  cross 
section  area  A  and  moment  of  inertia  I  may  vary  along  the 
length  but  do  not  alter  the  basic  assumption  that  plane 
cross  sections  remain  plane  after  loading. 

£.  All  cross  sections  share  an  initially  straight,  common  cen- 
troidal  axis;  the  x-axis  is  shown  in  Figure  1. 

d..  All  cross  sections  have  a  principal  axis  parallel  to  the 
y-axis.  Figure  1. 

£.  Deformations  due  to  shearing  strains  are  negligibly  small. 

f_.  Stress  concentrations  at  abrupt  changes  in  cross  sections, 
at  concentrated  loads  or  springs  and  at  supports  are  negli¬ 
gibly  smal 1 . 

Pi splacements 

6.  The  following  assumptions  regarding  displacements  of  the  system 
are  employed: 

a_.  Displacements  of  points  on  the  beam  centroidal  axis  are  com¬ 
pletely  defined  by  translations  u  (paral lei  to  the  x-axis), 
and  v  (parallel  to  the  y-axis),  and  rotation  9  about  an  axis 
perpendicular  to  the  x-y  plane. 


AM  displacements  are  small  and  do  not  alter  the  basic  geo¬ 
metry  of  the  system. 

£.  Rotations  8  are  small  enough  to  permit  the  approximations 
cos  8  =  1  and  sin  8  *  8 . 


died  Loads 


The  following  assumptions  regarding  external  applied  loads  are 
d: 

£.  Applied  loads  are  unaffected  either  in  magnitude,  direction 
or  point  of  application  by  the  deflections  of  the  system. 

b^.  Concentrated  applied  loads  are  assumed  to  be  applied  as  com¬ 
ponents  parallel  to  the  x  and  y  directions  or  as  couples 
acting  about  axes  perpendicular  to  the  x-y  plane. 

£.  Distributed  applied  loads  are  assumed  to  act  parallel  to  the 
x  and  y  axes. 


Fixed  Supports 


A  fixed  support  is  assumed  to  be  an  external  influence  which  re- 
n  a  specific  value  (either  zero  or  nonzero)  of  one  or  more  dis- 
nt  components  at  a  point  on  the  structure  regardless  of  other 


Concentrated  Linear  Spring  Supports 


Concentrated  linear  springs  produce  forces  on  the  structure 
re  directly  proportional  to  the  displacements  of  the  point  of 
int.  The  following  assumptions  are  employed: 
a.  Concentrated  linear  translational  springs  may  be  inclined 
at  an  angle  (a,  Figure  1)  to  the  axis  of  the  beam.  The  re¬ 


sisting  force  developed  in  the  spring  is  directly  propor¬ 


tional  to  the  translation  of  the  point  of  attachment  paral¬ 


lel  to  the  line  of  action  of  the  spring. 


J>.  A  concentrated  linear  rotational  spring  applies  a  resisting 
couple  proportional  to  the  rotation  of  the  point  of  attach¬ 
ment. 


Distributed  Linear  Spring  Supports 

10.  Distributed  linear  spring  supports  produce  distributed  resisting 
forces  directly  proportional  to  the  displacements  of  the  structure.  The 
following  assumptions  are  employed: 

£.  Distributed  linear  springs  resist  only  lateral  displacements 
v  or  axial  displacements  u. 

b^.  The  magnitude  of  the  distributed  resisting  force  depends 
only  on  the  lateral  displacement  v  or  axial  displacement  u 
at  any  point,  i.e.,  the  Winkler  hypothesis. 

Concentrated  Nonlinear  Spring  Supports 

11.  Concentrated  nonlinear  spring  supports  produce  forces  on  the 
structure  which  depend  on,  but  are  not  directly  proportional  to,  the  dis¬ 
placements  of  the  point  of  attachment.  The  following  assumptions  are 
employed: 

£.  Concentrated  nonlinear  translational  springs  may  be  inclined 
at  an  angle  (a.  Figure  1)  with  the  axis  of  the  beam.  The 
resisting  force  developed  in  the  spring  is  a  function  of  the 
translation  of  the  point  of  attachment  parallel  to  the  line 
of  action  of  the  spring. 

b^.  Nonlinear  rotational  springs  are  not  considered. 

Distributed  Nonlinear  Spring  Supports 

12.  Distributed  nonlinear  spring  supports  produce  distributed  re¬ 
sisting  forces  which  depend  on,  but  are  not  directly  proportional  to,  the 
displacements  of  the  structure.  The  following  assumptions  are  employed: 

£.  Distributed  nonlinear  springs  resist  only  lateral  displace¬ 
ments  v  or  axial  displacements  u. 


9 


b.  The  magnitudes  of  the  distributed  resisting  force  depend 
oniy  on  the  lateral  displacement  v  or  axial  displacement  u 
at  any  point,  i.e.,  the  Winkler  hypothesis. 

Characteristics  of  Nonlinear  Springs 

13*  The  force-deformation  relationship  for  a  nonlinear  spring  (ei ther 
concentrated  or  distributed)  is  assumed  to  be  provided  as  a  piecewise 
linear  curve,  as  shown  in  Figure  2.  For  any  deformation  & ,  parallel  to 
the  spring  line  of  action,  the  total  resisting  force  (concentrated  or  dis¬ 
tributed)  may  be  expressed  as 

Q  (or  q)  =  Q  (or  q  )  +  k&  (I) 

o  o 

Hence,  although  an  iterative  solution  is  required,  during  any  iteration  a 
nonlinear  spring  may  be  replaced  by  a  combination  of  a  fixed  load  (Qq  or 
qQ)  and  a  linear  spring  with  stiffness  k. 


Figure  2.  Characteristics  of  Nonlinear  Springs 


PART  III:  FINITE  ELEMENT  MODEL 


General 


l1*.  Solutions  for  the  general  beam-column  system  are  obtained  with 
the  finite  element  model  described  below. 

Nodes 


15*  Because  the  centroidal  axis  of  the  beam  is  initially  straight, 
the  location  of  a  node  is  completely  defined  by  its  x-coordinate.  For 
identification  in  the  following  paragraphs,  nodes  are  assumed  to  be  num¬ 
bered  sequentially  starting  with  node  I  at  the  left  end  of  the  beam,  as 
shown  in  Figure  3* 

16.  A  minimum  number  of  nodes  are  dictated  by  the  characteristics 
of  the  system: 

a_.  At  the  left  and  right  ends  of  the  beam. 

b^.  At  points  of  change  in  beam  properties  (i.e.,  E,  I,  and/or 
A). 

£.  At  the  point  of  application  of  concentrated  loads. 

<£.  At  fixed  supports. 

£.  At  the  point  of  application  of  concentrated  spring  supports. 

f_.  At  the  beginning  and  end  of  distributed  (either  axial  or 
lateral)  loads. 

£.  At  the  beginning  and  end  of  distributed  spring  (either  axial 
or  lateral)  supports. 

17.  Displacements,  internal  forces,  and  other  effects  are  determined 
only  at  the  nodes.  To  provide  more  detailed  description  of  the  variation 
of  effects  along  the  length,  nodes  may  be  defined  at  any  other  location  on 
the  beam.  In  addition,  the  accuracy  of  a  solution  may  be  affected  by  the 
number  of  nodes.  This  aspect  will  be  discussed  subsequently. 

Variations  in  System  Properties 

18.  The  following  assumptions  are  employed  for  system  properties 
which  vary  along  the  length  of  the  structure. 


II 


a.  Beam  cross  section  properties  E,  A,  and  I  are  assumed  to  be 
constant  between  adjacent  nodes.  In  effect,  a  tapered  beam 
is  replaced  by  a  piecewise  prismatic  (stepped)  member.  As 
discussed  later  this  replacement  is  accomplished  automatic¬ 
al  ly . 

b^  Distributed  loads  (either  applied  or  resulting  from  the 
fixed  load  portion  of  distributed  nonlinear  springs)  are 
assumed  to  vary  linearly  between  adjacent  nodes. 

£.  Distributed  spring  stiffnesses  (either  distributed  linear 
springs  or  the  linear  stiffness  component  of  a  distributed 
nonlinear  spring)  are  assumed  to  vary  linearly  between  adja 
cent  nodes. 


Elements 


19.  An  element  is  defined  as  the  portion  of  the  beam  between  adja¬ 
cent  nodes.  For  reference  in  the  following  paragraphs,  an  element  is 
identified  by  the  node  number  at  its  left  end,  as  shown  in  Figure  3. 

20.  Under  the  assumptions  stated  above,  each  element  is  character¬ 
ized  by  modulus  of  elasticity  E,  cross  section  properties  A  and  I,  length 
h,  and  is  subjected  to  distributed  loads  and  springs,  as  shown  in  Figure 
A. 

21.  Interaction  of  internal  axial  stress  resultants  with  lateral  de¬ 
flections  will  affect  the  bending  response  of  the  element.  When  this 
effect  is  included,  the  bending  resistance  of  the  element  is  developed 
using  a  constant  average  axial  stress  resultant  given  by 

P*lf{ui+1  '  ui}  (2) 

22.  The  relationship  between  element  end  forces  and  end  nodal  dis¬ 
placements,  obtained  by  procedures  given  in  References  (1),  (2),  and  (3) 
is 

f5  “  [s'E  +  S^J  u'  +  f'e  (3) 
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where 


Distributed 

Forces 


1  \W 

W,  A^|kx>1+1 

CD 

- - - - - 1 

: _ i  1  ky.i*' 

Distributed 

Springs 

Figure  k.  Typical  Element 


fl  -  f!A,iT  -  tf'  .  f1  .  «!  i  f!  ...  f1  ...  ,iT 

“i  ~j  +  )J  s,i  y j i  i  |  x,i+l  y,i+l  i+l 

=  (6x1)  vector  of  element  end  forces; 

Sg  =  (6x6)  axial  and  beam  bending  stiffness  matrix  including  effects 


of  axial  stress  resultant  on  bending; 

(6x6)  stiffness  matrix  representing  effects  of  distributed 
springs; 


u'  =  [U.  U.  .]T  =  [u.  v.  e.  ;  u. 
i  ~i+l  i  i  r  J  i 


+1  Vi+1  9i+l^ 


=  (6x1)  vector  of  end  nodal  displacements;  and 


t  1  =  if* .  f:  -j 


[f‘  .  f'  .  M 1  .  J  f' 

AV  I  AW  I  A  I  '  C 


M 


J 


]T  _ 

•e,i  ~e,i+l  1  exfi  ey,i  e,i  J  ex, i+l  ey,i+l  e,i+l 

=  (6x1)  vector  of  fixed  end  forces  due  to  distributed  loads. 

i  i  1 

23-  Coefficients  for  matrices  §g,  S^,  and  f  are  given  in  Figures 
5,  6,  and  7  end  Table  1. 


Node  Equi 1 ibrium 


24.  A  typical  node,  shown  in  Figure  8,  is  subjected  to: 

a.  Concentrated  forces  F  .,  F  .,  and  C.,  where  F  .  and  F 
—  x, i  y, i  i  x, i  y , t 

include  both  applied  external  concentrated  loads  and  the 
fixed  force  component  of  all  concentrated  nonlinear  springs 
attached  to  the  node. 

b^.  One  or  more  concentrated  linear  translation  springs  which 
include  the  linear  spring  components  of  concentrated  non¬ 
linear  springs  attached  to  the  node. 

£.  One  or  more  linear  rotation  springs. 

d^.  Element  end  forces  from  elements  on  either  side  of  the  node. 

25.  Equilibrium  of  element  i  is  expressed  by 


k  .cos  a  k  . sinacosa 
c,t  c,i 


k  .sinacosa  k  .sin  a 
c,i  c,i 


0 


0 


\ 


0 


iVi 


i 


or 


f!*1  +  f!  + 

-  I  ■*’  I 


F. 


(4) 
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Figure  5-  Element  Axial  and  Beam  Bending  Stiffness  Matrix 
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■  constants  given  in 
Table  1 


Figure  7.  Element  Fixed  End  Forces  Due  to  Distributed  Loads 
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where  the  summation  Z  extends  over  all  concentrated  springs  attached  to 
to  node  i . 

26.  Substitution  for  end  force  vectors  f!  and  f!  in  terms  of  ele¬ 
ment  stiffness  matrices  using  Equation  (3)  and  Figures  5,  6,  and  7  results 
in  the  following  governing  equation  which  must  be  satisfied  at  every  node. 


re'-1  +  e 1 ” 1 

~  E ,  i ,  i  - 1  'K,i,i-lJ  - i-1 


.  i  - 1 


,]  u: 


+  [s'"!  .  +  s'"!  .  +  s'  .  .  +  s'  .  .  +  ZK  .]  U. 

^  E , »  ,  i  *”C,i  'i 


+  [s'  .  .  .  +  s'  .  .  U.  .  *  F.  +  f'"!  +  f'  . 

-E,i,i+1  -K, i,i+l  'i+l  '•  'e,i  'e,i 


(5) 


27-  Evaluation  of  Equation  (5)  at  every  node  in  the  model  results  in 
a  system  of  simultaneous  equations  (3n  equations  for  a  model  with  n  nodes) 
which,  after  revision  to  reflect  fixed  supports,  may  be  solved  for  nodal 
displacements.  When  nodal  displacements  are  known,  internal  forces, 
spring  forces,  and  reactions  may  be  determined  for  the  system. 


I  terat ion 


28.  A  single  solution  of  the  simultaneous  equations  produces  final 
results  for  systems  possessing  the  following  characteristics: 

£.  Only  li  .^ar  spring  supports  are  present;  and 
b^.  The  effect  of  axial  stress  resultant  on  element  bending 
stiffness  is  excluded;  or 

£.  The  effect  of  axial  stress  resultant  on  element  bending 

stiffness  is  included  and  axial  distributed  springs  are  ab¬ 
sent  and  all  concentrated  translational  springs  are  perpen¬ 
dicular  to  the  beam  axis. 

29.  All  other  systems  require  an  iterative  solution  in  order  to 
evaluate  axial  stress  resultants  and/or  fixed  force-linear  spring  compo¬ 
nents  for  nonlinear  springs.  The  iterative  solution  process  is  initiated 
by  evaluating  axial  stress  resultants  and  nonlinear  spring  characteris¬ 
tics  for  zero  displacements.  The  simultaneous  equations  are  evaluated 
and  solved  for  a  new  estimate  of  system  displacements.  On  each  iteration 
the  solution  of  the  preceding  iteration  is  used  for  evaluation  of  system 
properties.  This  process  is  continued  until  the  results  of  two  successive 
iterations  differ  by  an  acceptably  small  amount. 
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Effect  of  Node  Spacing  on  Solution 

.  The  finite  element  procedure  described  above  produces  "exact" 
ons  for  systems  possessing  the  following  characteristics: 

a.  The  beam  is  piecewise  prismatic. 

b^  Only  concentrated  external  loads  and  linearly  varying 

lateral  (y)  loads  are  present.  Or  if  the  effect  of  axial 
stress  resultant  on  element  bending  stiffness  is  excluded, 
distributed  axial  (x)  loads  vary  linearly. 

£.  Only  concentrated  springs  are  present  and,  if  the  effect  of 
axial  stress  resultant  on  element  bending  stiffness  is  in¬ 
cluded,  all  concentrated  translation  springs  are  perpendicu¬ 
lar  to  the  beam  axis.  (Note:  The  accuracy  of  a  solution 
including  the  effect  of  axial  stress  resultant  in  the  pre¬ 
sence  of  inclined  translational  springs  may  be  affected  by 
the  iteration  convergence  tolerance  but  not  by  node  spac¬ 
ing.) 

d.  Adjacent  elements  are  approximately  equal  in  length.  Adja¬ 
cent  elements  having  drastically  different  lengths  may 
result  in  significant  round-off  errors  in  the  solution 
of  simultaneous  equations. 

In  all  other  systems  the  number  of  nodes  (and  elements)  used  in 
nite  element  model  affects  the  accuracy  of  the  solution.  In  gener- 
the  number  of  nodes  (and  elements)  is  increased,  the  solution 
to  converge  to  the  "exact"  solution.  There  is  no  "rule-of-thumb" 
i/ide  the  necessary  number  of  nodes  for  acceptable  results.  It  may 
issary  to  perform  several  solutions  for  various  numbers  of  nodes 
Je  spacings  to  ensure  that  an  adequate  solution  has  been  obtained. 
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PART  IV:  COMPUTER  PROGRAM 


General 


32.  A  computer  program  implementing  the  analytical  process  described 
above  has  been  written  in  the  FORTRAN  programming  lanquage.  With  minor 
revisions  the  program  will  be  operational  on  a  computer  employing  a  60  bit 
('15  decimal  digit)  word  length.  For  systems  using  fewer  than  15  decimal 
digits  it  may  be  necessary  to  perform  some  arithmetic  operations  in  double 
preci s ion . 

Input  Data 

33-  The  program  has  been  written  to  generate  automatically  intermedi¬ 
ate  data  values  from  a  minimum  of  user  input  data.  Details  of  required 
user  input  data  are  presented  in  the  Input  Guide,  Appendix  A.  The  pro¬ 
cesses  and  assumptions  used  in  converting  user  input  data  to  intermediate 
data  values  are  described  below. 

Global  Coordinate  System 


3A.  A  horizontal  beam  centroidal  axis  (the  x-axis,  positive  to  the 
right)  is  assumed.  The  global  y  axis  is  positive  upward.  The  origin  of 
the  x-y  coordinate  system  is  arbitrary.  A  local  coordinate  system  asso¬ 
ciated  with  nonlinear  springs  is  described  later. 

Displacement  and  Load  Sign  Conventions 

35-  Positive  directions  for  displacements  and  loads  are  as  follows: 
£.  Nodal  displacements  u  and  v  are  positive  if  the  node  trans¬ 
lates  in  the  positive  x  and  y  directions,  respectively. 
b_.  Nodal  rotation  is  positive  if  the  node  rotates  counterclock¬ 
wise. 

c.  External  loads  (concentrated  F  and  F  :  or  distributed  q 
—  x  y  ^x 

and  q^)  are  positive  if  they  act  in  the  positive  x  and  y 
directions,  respectively. 


d_.  A  concentrated  applied  couple  (C)  is  positive  if  it  tends  to 
rotate  the  node  counterclockwise. 

Concentrated  or  distributed  linear  spring  stiffnesses  are  al¬ 
ways  positive. 


Nonlinear  Spring  Conventions 


36.  As  stated  previously,  the  characteristics  of  nonlinear  springs 
are  assumed  to  be  described  by  a  piecewise  linear  curve  giving  spring  re¬ 
sisting  force  as  a  function  of  spring  deformation.  The  sign  conventions 
assumed  for  each  of  the  three  types  of  nonlinear  springs  permitted  by  the 
program  are  as  follows. 


Concentrated  Nonlinear  Springs 


37.  Concentrated  nonlinear  springs  are  assumed  to  resist  nodal  trans¬ 
lation  displacement  components  parallel  to  the  spring  line  of  action. 
Positive  resisting  forces  and  spring  deformations  for  concentrated  springs 
are  defined  for  a  local  x1  coordinate,  as  shown  in  Figure  9-  The  origin 
of  the  local  x'  coordinate  is  at  the  point  of  attachment  (node);  the  posi¬ 
tive  x1  direction  extends  along  the  spring  line  of  action  away  from  the 
point  of  attachment.  The  inclination  of  the  spring  is  given  by  the  angle 
a  between  the  global  x-axis  and  the  local  x'  axis  (positive  counterclock¬ 
wise)  . 

38.  Spring  resisting  forces  are  positive  if  they  act  in  the  positive 
local  x'  dire:  on.  To  illustrate,  consider  a  spring  in  an  initial  state 
o!  compression  at  zero  nodal  displacement.  The  spring  therefore  would 
prod<c.  force  act  i  ng  on  the  node  in  the  negative  x'  direction  (i.e., 
negative  force),  -Fj ,  at  zero  spring  deformation,  A^  =  0,  as  illustrated 
by  point  (I)  in  Figure  9b.  If  the  point  of  attachment  undergoes  nodal 
displacements  u  and  v  to  produce  a  positive  deformation  +A^,  Figure  9a, 
the  compression  (i.e.,  negative  spring  force),  -F^,  would  increase  as 
illustrated  by  point  (2)  in  Figure  9b.  Further  positive  deformation  to 
Aj  would  result  in  a  negative  resisting  force,  -F^,  as  shown  by  point  (3). 
If  the  nodal  displacements  produce  deformation  in  the  negative  x'  direc¬ 
tion,  e.g.,  -A^,  the  compressive  (negative)  force  in  the  spring  would 
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a.  Concentrated  Nonlinear  Spring 


Spring  Resistance 
in  x1  Direction 


Spring  Deformation 
in  x1  Direction 


b.  Resisting  Force-Deformat  ion  Curve 
Figure  9.  Characteristics  of  Concentrated  Nonlinear  Spring 
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reduce  to,  -F^,  point  (4).  At  some  negative  deformation,  -A,.,  the  ini¬ 
tial  compression  might  be  totally  relieved,  i.e.,  F,.  =  0,  as  shown  by 
point  (5).  Further  negative  deformation,  -A^,  might  produce  tension  in 
the  spring  (+F^),  as  illustrated  by  point  (6).  The  coordinate  values, 
A,F,  of  each  point  on  the  deformation-resistance  curve  are  provided  as 
input.  A  linear  variation  between  input  deformation-resistance  coordi¬ 
nates  is  assumed. 

Distributed  Nonlinear  Springs 

39.  Distributed  nonlinear  springs  produce  distributed  forces  propor¬ 
tional  to  the  translation  displacements  of  the  nodes.  The  deformation  of 
the  spring  is  equal  to  and  has  the  same  sign  as  the  nodal  displacement. 
The  distributed  force  produced  on  the  beam  by  a  distributed  nonlinear 
spring  is  positive  if  the  force  acts  in  the  positive  global  x  or  y  coor¬ 
dinate  directions.  Distributed  spring  characteristics  are  assumed  to  be 
described  by  pairs  resisting  force-deformation  coordinate  values  similar 
to  those  for  concentrated  springs. 

Data  Generation 


40.  Concentrated  applied  loads,  concentrated  springs  and  fixed  sup¬ 
ports  occur  at  isolated  points  on  the  structure.  All  other  data  (beam 
properties  E,  A,  and  I,  and  distributed  loads  and  springs)  are  distribut¬ 
ed  over  a  range  of  x-coord i nates .  The  user  is  required  to  provide  the 
locations  (x-coordinates)  at  which  concentrated  effects  occur  and  the 
terminal  values  (beginning  and  end)  of  distributed  quantities.  These 
define  the  "dictated"  nodes  in  the  model.  The  program  also  provides  for 
defining  nodes  intermediate  to  the  dictated  nodes.  The  following  para¬ 
graphs  describe  the  procedures  and  assumptions  used  to  obtain  system 
characteristics  at  the  intermediate  nodes. 


Beam  Cross  Section  Properties 

41.  Beam  cross  section  properties  A  and  I  are  provided  at  the  ter¬ 
minals  of  a  distribution.  Figure  10  illustrates  for  cross  section  area 
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the  steps  used  to  convert  input  data  to  element  properties.  The  user  has 

specified  values  of  cross  section  area  A.  and  A.+i<  at  x.  and  x.+i<,  respec 

tively,  resulting  in  "dictated"  nodes  i  and  i+4.  Through  node  spacing 

data  (see  below)  intermediate  nodes  have  been  defined  at  i+1,  i+2,  and 

i+3.  The  initial  assumption  is  made  that  cross  section  area  varies 

linearly  from  A.  at  x.  to  A.+^  at  x.+lj,  as  shown  in  Figure  10b,  to  result 

in  values  of  cross  section  areas  A...,  A.,.,  and  A.,,  at  x...,  x..„,  and 

i+l  i+2  i+3  i+l  i+2 

x.+j,  respectively.  The  distribution  is  further  converted  to  the  piece- 
wise  prismatic  representation  shown  in  Figure  10c  by  assigning  to  each 
element  covered  by  the  distribution  a  constant  area  equal  to  the  average 
of  the  areas  at  each  end  of  the  element,  e.g.,  A,+  =  (A.+j  +  h.^+^/2. 
Note  that  E  is  required  to  be  constant  over  the  distribution. 

42.  Beam  cross  section  properties  data  are  assumed  to  define  the  ex¬ 
treme  (left  and  right)  ends  of  the  beam-column  system  to  be  analyzed. 
Nonzero  values  of  section  properties  E,  A,  and  I  must  be  supplied  for 
every  point  on  the  structure  between  these  extremes.  Any  other  distri¬ 
buted  or  concentrated  data  which  fall  beyond  these  extremes  are  ignored 
by  the  program. 


Node  Spacing  Data 

43.  Node  spacing  data  provide  for  defining  nodes  intermediate  to  the 
"dictated"  nodes.  The  user  specifies  a  maximum  node  spacing, H  ,  to  be 

(H3X 

used  in  a  range  of  x-coord i nates .  For  example,  consider  the  range  of  x- 
coordinates  between  "dictated"  nodes  x.  and  x.+i>  in  Figure  10a.  The 
range  of  coordinates  x.  to  x.+i{  will  be  subdivided  into  elements  of  equal 
length  as  follows.  The  number  of  elements  is  equal  to  the  largest  integer 
given  by 

<»l*  -  <6) 

and  the  length  of  elements  in  this  region  will  be 

h  -  (x.+Jj  -  x.)/n  (7) 

If  is  greater  than  or  equal  to  (x.+i>  _  x.),  then  the  entire  region 

will  be  treated  as  a  single  element. 


44.  If  a  "dictated"  node  falls  between  the  limits  of  a  range  of  node 
spacing  data,  intervals  to  the  left  and  right  of  the  dictated  node  are 
treated  as  separate  ranges  for  intermediate  node  spacing. 

Distributed  Loads 


45-  Distributed  loads  are  assumed  to  vary  linearly  between  the  values 
at  the  beginning  and  end  of  the  distribution  provided  as  input. 

Distributed  Linear  Springs 

46.  The  program  provides  for  variations  in  linear  spring  stiffness 
coefficients  over  an  interval  to  given  by  the  general  expression 

k  =  a  +  bzC  (8) 

where 

k  =  stiffness  coefficient  (force/length/displacement)  of  distributed 
x  or  y  spring  at  any  point  x^  to  Xjl 
a  =  stiffness  at  always  positive; 

Q 

b  =  constant  with  units  such  that  bz  has  units  of  stiffness;  b  may 
be  either  positive  or  negative;  however,  only  positive  values  of 
k  are  permitted; 

z  »>  distance  measured  from  and 
c  =  dimensionless  constant. 

47.  For  the  input  values  x.  ,  x.+^,  a,  b,  and  c,  the  program  calcu¬ 
lates  at  each  intermediate  node  the  distributed  spring  stiffness  from  the 
general  expression  given  in  Equation  (8),  as  illustrated  in  Figure  11. 

For  analysis,  the  nonlinear  distribution  given  by  the  general  Equation  (9) 
is  replaced  by  the  piecewise  linear  variation  shown  by  the  dashed  lines  in 
F i gure  1  I . 


Distributed  Nonlinear  Springs 


48.  Distributed  nonlinear  springs  are  described  by 
deformation  curves  for  each  end  of  the  distribution,  as 
Figure  12.  As  stated  previously,  a  nonlinear  spring  is 


resisting  force- 
illustrated  in 
represented  by  a 
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combination  of  a  fixed  force  and  a  linear  spring.  The  components  of  the 

nonlinear  spring  at  intermediate  nodes  are  obtained  as  follows.  Consider 

node  i+l  in  Figure  12.  The  spring  deformation  at  node  i+l,A.+j  (either 

II.,,  or  v.,.)  is  used  with  the  input  curve  at  node  i  to  obtain  q  .  and 
l+l  l+l  0,4 

k. ;  and  with  the  curve  at  node  i+3  to  obtain  q  and  k  .  The  values  at 
l  o,r  r 

node  l+l  are  then  calculated  from  linear  interpolation  as 


Ho, i+l 


l%,i  (xi+3  '  Xi+1)  +  %,r  (xi  +  l  ‘  xi)1/(xi+3  ‘  xi 


'9) 


and 


i  +  l 


[k*  (xi+3  ‘  xi+l5  +  kr  (xi+l  '  xi)1/(xi+3 


xi) 


10) 


Output  Data 

k3.  Output  data  are  provided  in  three  parts.  Output  may  be  directed 
to  the  user  terminal,  to  a  data  file,  to  both  or  some  parts  of  the  output 
may  be  omitted  entirely. 


Echoprint  of  Input  Data 

50.  This  section  contains  a  tabular  listing  of  all  input  data  cur¬ 
rently  available  to  the  program.  This  section  is  output  at  the  option  of 
the  user. 


Summary  of  Results 


51.  This  section  contains:  maximum  positive  and  maximum  negative 
values  of  displacements  and  internal  forces  and  the  x-coordinates  at  which 
the  maxima  occur;  reactions  at  fixed  supports;  and  forces  in  concentrated 
linear  and/or  nonlinear  springs.  This  section  is  always  output  to  either 
the  user  terminal,  to  an  output  file,  or  both. 

Complete  Results 

52.  This  section  contains  a  tabulation  of  displacements,  internal 
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stress  resultants,  and  forces  in  distributed  springs  at  nodes  contained 
within  a  range  of  x-coordinates  specified  by  the  user.  This  section  may 
be  omitted  entirely. 


Sign  Conventions  for  Output 


53.  Sign  conventions  used  for  output  data  are  presented  in  Table  2. 


Table  2.  Output  Sign  Conventions 


Quant i ty 


Sign  Convention 


Translations  u  and  v 


Rotation  6 

Axial  Internal  Stress 
Resultant 

Shear  Stress  Resultant 


Bending  Moment 

X,Y  Reactions  at  Fixed 
Supports 

Moment  Reactions  at 
Fixed  Supports 

Forces  in  Concentrated 
Linear  Translation 
Springs 

Moment  in  Linear  Rota¬ 
tion  Spring 

Forces  in  Concentrated 
Nonlinear  Translation 
Springs 

Forces  in  Distributed 
Linear  Springs 

Forces  in  Distributed 
Nonlinear  Springs 


Positive  if  node  moves  in  positive  x  or  y  coor¬ 
dinate  direction,  respectively 

Positive  if  node  rotates  counterclockwise 

Positive  if  tension 


Positive  if  shear  force  tends  to  move  left  end 
of  segment  of  beam  in  positive  y  direction 
(i.e.,  up  on  left,  down  on  right  end  of  segment) 

Positive  if  produces  compression  in  plus  y 
face  (top)  of  beam 

Positive  if  act  in  positive  x  and  y  coordinate 
directions,  respectively 

Positive  if  act  counterclockwise  on  beam 
Positive  if  tension 


Positive  if  acts  counterclockwise  on  beam 

Sign  determined  according  to  input  resisting 
force-deformation  curve  coordinates 

Positive  if  tension 

Sign  determined  according  to  input  resting 
force-deformation  curve  coordinates 


PART  V:  EXAMPLE  SOLUTIONS 


General 


5*t.  Program  output  and  supporting  information  demonstrating  the  use 
of  the  program  are  contained  in  Appendix  B.  Principal  features  of  each 
solution  are  discussed  below. 

55.  Interactive  control  of  the  program  is  shown  with  the  output  for 
Example  I. 


Example  1:  Fixed  End  Beam 

56.  The  fixed  end  beam  shown  in  Figure  B1  was  analyzed  for  a  succes¬ 
sion  of  loading  conditions  described  in  the  input  data  file  shown  on  page 
B3. 

a_.  The  basic  problem  is  described  in  lines  1000  through  1120  of 
the  file.  During  execution,  effects  of  axial  load  on  bending 
were  excluded  for  the  solution  of  this  part. 

b^.  On  the  first  rerun.  Example  1A,  only  the  problem  heading  was 
changed  to  that  shown  on  lines  11 1*0  through  11 60  of  the  data 
file.  All  other  data  remain  unchanged.  During  execution 
the  effects  on  bending  due  to  the  40kip*  axial  load  were  in¬ 
cluded  in  the  solution. 

£.  On  subsequent  reruns,  lines  1170  through  1340,  the  axial 

load  was  increased  by  adding  concentrated  load  data  to  each 
•preceding  data  set. 

57.  Echoprints  of  input  data  were  requested  only  for  the  first  and 
last  in  the  succession  of  problems,  pages  B5  and  BI2,  respectively. 

58.  The  summary  of  results  was  output  for  each  problem.  Because  the 
model  used  to  represent  the  beam  contains  only  three  nodes  and  two  ele¬ 
ments,  no  additional  information  would  have  been  provided  by  a  complete 
tabulation. 

59.  The  deflections  and  forces  presented  are  "exact"  solutions  for 
each  loading.  The  effect  of  axial  load  on  bending  is  to  increase  the  de¬ 
flections  and  bending  moments.  As  the  total  axial  load  approaches  the 
theoretical  buckling  load  (3289.87  kips),  deflections  and  moments  attain 

*  A  table  of  factors  for  converting  inch-pound  units  of  measurement  to 
metric  (SI)  units  is  presented  on  page  4. 
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extremely  large  values.  However, for  an  axial  load  below  buckling,  deflec¬ 
tions  and  moments  are  consistent  with  the  direction  of  the  applied  later¬ 
al  load.  For  an  axial  load  slightly  greater  than  the  buckling  load,  de¬ 
flections  and  moments  undergo  a  reversal  in  sign  (i.e.,  are  inconsistent 
with  lateral  load  direction)  indicating  that  buckling  has  occurred. 

Example  2:  Beam  on  Uniform  Elastic  Foundation 

60.  Accuracy  of  results  produced  by  the  finite  element  model  used  to 
represent  the  system  shown  in  F igure  82  depends  on  node  spacing.  Solutions 
were  obtained  for  two  node  spacings  as  indicated  in  the  data  file  shown  on 
page  BI5. 

61.  The  echoprint  of  input  data,  summary  of  results,  and  the  (trivi¬ 
al)  complete  results  table  for  the  first  case  (lines  1000  through  11  AO  of 
the  input  file)  are  shown  on  pages  BI6  and  BI7.  Maximum  lateral  deflec¬ 
tion  and  bending  moment  for  the  model  wi th  three  nodes  and  two  elements 
differ  only  slightly  from  the  "exact"  sotution  (0. 19097  in.  and -1 .52025E+5, 
respectively) . 

62.  Increasing  the  model  to  eleven  nodes  and  ten  elements  produces 
results,  pages  Bl  8  and  Bl  9,  whi ch  are  exact  to  four  significant  figures. 

63.  It  should  be  emphasized  that  node  spacing  may  have  a  more  signi¬ 
ficant  effect  on  the  accuracy  of  the  solution  than  is  indicated  by  this 
problem.  This  is  demonstrated  in  Example  3- 

Example  }:  Pile  Head  Stiffness  Matrix 

6A.  The  program  was  used  to  generate  the  pile  head  stiffness  matrix 
for  the  system  shown  in  Figure  B3.  By  definition,  coefficients  of  the 
pile  head  stiffness  matrix  are  the  reactions  generated  by  imposing  unit 
values  of  deflection  components  at  the  pile  head. 

65.  The  problem  was  run  using  the  input  data  file  shown  on  page  B2I. 
Because  axial  load  effects  on  bending  were  excluded,  axial  and  lateral 
coefficients  for  the  stiffness  matrix  are  obtained  simultaneously  by  im¬ 
posing  unit  axial  (x)  and  lateral  (y)  deflections  at  the  pile  head  with 
rotation  at  that  point  equal  to  zero.  Two  solutions  are  presented  for 
this  case:  for  nodes  spaced  at  5  ft,  and  nodes  spaced  at  1  ft.  Moment 
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coefficients  for  the  pile  head  stiffness  matrix  were  obtained  for  a  I  ft 
node  spacing  by  specifying  zero  axial  and  lateral  displacements  and  a 
unit  rotation  at  the  pile  head. 

66.  Results  for  each  case  are  given  on  pages  022  through  B27.  Values 
of  stiffness  coefficients  obtained  for  the  5  ft  node  spacing  are  signifi¬ 
cantly  in  error.  Decreasing  the  node  spacing  to  1  ft  produced  the  pile 
head  stiffness  matrix  shown  below.  Values  shown  in  parentheses  were  ob¬ 
tained  by  another  approximate  method  reported  in  Reference  (4). 
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67-  Reducing  the  node  spacing  to  0.5  ft  had  little  effect  on  the  re¬ 
sults. 


Example  4:  Anchored  Retaininc 


68.  A  one  foot  strip  of  the  multiple  anchored  retaining  wall  describ¬ 
ed  in  Figure  B4  was  analyzed  as  a  beam-column  in  which  the  soil  and  anchors 
were  represented  by  nonlinear  springs. 

69.  Procedures  described  in  Reference  (5)  were  used  to  obtain  the 
force-displacement  curves  shown  in  Figure  B5  for  the  soil  resisting  latei — 
al  (y)  displacements  of  the  system.  To  the  writer's  knowledge  no  proced¬ 
ure  has  been  established  for  representing  wall-soil  friction  interaction. 
The  force-d i sp lacement  curves  shown  in  Figure  B6  were  arbitrarily  selected 
for  this  effect. 

70.  The  anchors  were  represented  by  concentrated  nonlinear  springs 
with  characteristics  shown  in  Figure  B7- 

71.  The  input  data  file  used  for  solution  with  the  program  is  shown 
on  page  B32.  Output  from  the  program  is  given  on  pages  B33  through  B38. 
Anchor  forces  are  given  by  section  II.C  of  the  output.  Soil  pressures 
are  given  by  the  tabulation  of  forces  in  distributed  nonlinear  springs, 
section  (II.C  of  the  output. 
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APPENDIX  A:  GUIDE  FOR  DATA  INPUT 


Source  of  Input 


1.  Input  data  may  be  supplied  from  a  predefined  data  file  or  from 
the  user  terminal  during  execution.  If  data  are  supplied  from  the  user 
terminal,  prompting  messages  are  printed  to  indicate  the  amount  and 
character  of  data  to  be  entered. 

Data  Ed i t i ng 

2.  When  all  data  for  a  problem  have  been  entered,  the  user  is 
offered  the  opportunity  to  review  an  echoprint  of  the  currently  avail¬ 
able  input  data  and  to  revise  any  or  all  sections  of  the  input  data  be¬ 
fore  execution  is  attempted.  When  editing  during  execution,  each  sec¬ 
tion  must  be  entered  in  its  entirety. 

Input  Data  File  Generation 

3.  After  data  have  been  entered  from  the  terminal,  either  initial¬ 
ly  or  after  editing,  the  user  may  direct  the  program  to  write  the  input 
data  to  a  permanent  file  in  input  data  file  format. 

Data  Format 


1*.  All  input  data  (whether  supplied  from  the  user  terminal  or  from 
a  file)  are  read  in  free  field  format: 

£.  Data  items  must  be  separated  by  one  or  more  blanks  (comma 
separators  are  not  permitted). 

1).  Integer  numbers  must  be  of  form  NNNN. 

c_.  Real  numbers  may  be  of  form 

±xxxx,  ±xx.xx,  or  ±xx.xxE+ee 

d..  User  responses  to  all  requests  for  control  by  the  program 
for  alphanumeric  input  may  be  abbreviated  by  the  first  let¬ 
ter  of  the  indicated  word  response,  e.g.  , 

ENTER  'YES'  OR  'NO'— respond  Y  or  N 
ENTER  'CONTINUE'  OR  ' END ' --respond  C  or  E 
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Carriage  return  responses  aione  will  cause  abnormal  termina 
tion  of  the  program. 


Sections  of  Input 

5-  Input  data  are  divided  into  the  following  sections: 

£.  Heading. 

b_.  Beam  Cross  Section  Properties. 

c_.  Node  Spacing. 

d..  Applied  Loads. 

e_.  Fixed  Supports. 

f_.  Linear  Spring  Supports. 

Nonlinear  Spring  Supports. 
h_.  Termination. 

6.  When  data  are  entered  from  the  terminal,  data  sections  a.  through 
g_ may  be  input  in  any  order.  When  data  are  entered  from  a  predefined 
data  file,  the  heading  section  £  must  be  entered  first;  other  sections  b^ 
through  £  may  be  entered  in  any  ord*-r. 

7.  When  data  are  entered  from  the  terminal,  the  user  is  cued  for 
Termination.  Termination  for  a  data  file  is  discussed  later. 

Minimum  Required  Data 

8.  Data  sections  £,  b_,  and  c_  are  always  required.  At  least  one  of 
sections  e_,  f,  and  £  is  required.  It  is  the  responsibility  of  the  user 
to  ensure  that  sufficient  supports  (either  fixed  or  spring  supports)  are 
provided  to  inhibit  all  rigid  body  motions  of  the  system,  i.e.,  to  ensure 
a  stable  structure. 


Uni  ts 


9-  The  program  recognizes  the  following  units: 

Jnches  Feet  Pounds  Kjps 

Default  is  to  inches  and  pounds.  Angular  units  are  either  degrees  or 
radians  as  explained  below. 

10.  Each  data  section  may  be  entered  with  any  combination  of  units 
for  length  and  force  as  specified  by  the  user. 
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Predefined  Data  File 


11.  In  addition  to  the  general  format  requirements  given  in  para¬ 
graph  4,  above,  the  following  pertain  to  a  predefined  data  file  and  to 
the  input  data  description  which  follows: 

a_.  Each  line  must  commence  with  a  nonzero,  positive  line  num¬ 
ber,  denoted  LN  below. 

b_.  A  line  of  input  may  require  both  alphanumeric  and  numeric 
data  items.  Alphanumeric  data  items  are  enclosed  in  sinal^ 
quotes  in  the  following  paragraphs. 

£.  A  line  of  input  may  require  a  keyword.  The  acceptable 
abbreviation  for  the  keyword  is  indicated  by  underlined 
capital  letters.  E.g.,  the  acceptable  abbreviation  for 
the  keyword  1 L I  Near 1  is  'LIN*. 

<d.  Lower  case  words  in  s;ngle  quotes  indicate  a  choice  of  key¬ 
words  defined  following. 

£.  I  terns  designated  by  upper  case  letters  and  numbers  without 
quotes  indicate  numeric  data  values.  Numeric  data  values 
are  either  real  or  integer  according  to  standard  FORTRAN 
variable  naming  conventions. 

f_.  Data  items  enclosed  in  brackets  [  ]  may  not  be  required. 

Data  items  enclosed  in  braces  {  }  indicate  special  note 

f o 1 1 ows . 

£.  Input  data  are  divided  into  the  sections  discussed  in  para¬ 
graph  5,  above.  Except  for  the  heading,  each  section  con¬ 

sists  of  a  header  line  and  one  or  more  data  lines.  The 
header  line  serves  the  multiple  porposes  of:  indicating 
the  end  of  the  preceding  section;  identifying  the  data  sec¬ 
tion  to  follow;  indicating  whether  the  section  is  'New'  or 
to  be  'Added'  to  the  same  section  of  a  preceding  problem 
(see  below);  and  indicating  the  units  associated  with  the 
section. 

h_.  A  data  file  may  contain  input  data  for  several  problems  to 
be  run  in  succession.  On  the  first  problem  defined  in  the 
file,  all  data  sections  must  be  'New'.  In  subsequent  prob¬ 
lem  descriptions,  'New'  indicates  the  section  is  to  replace 
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that  section  in  the  preceding  problem;  'Add'  indicates  the 
data  are  to  be  appended  to  that  section  in  the  preceding 
problem.  A  data  section  may  be  deleted  from  the  preceding 
problem  description  by  indicating  'New'  on  the  section  head¬ 
er  line  and  immediately  following  with  another  section  head¬ 
er. 

j_.  Each  section  header  line  contains  the  information  [{'New1  or 
'Add'}]  and  [{'units'}].  Either  or  both  of  these  items  may 
be  omi tted . 

_i_(  1 )  If  [{'New'  or  'Add'}]  is  omitted,  'New'  is  assumed. 

Omit  this  item  on  the  first  problem  description  in  the 
file. 

j_(2)  If  'Add'  is  indicated,  the  units  for  appended  data 

must  be  the  same  as  for  that  section  in  the  preceding 
problem. 

_[_( 3)  [{'units'}]  must  be  'inches'  or  'Feet'  and  'Pounds'  or 
•J<i ps ' .  If  data  are  to  be  appended  to  the  preceding 
problem,  omit  this  item.  Data  to  be  appended  must  be 
in  the  same  units  prescribed  for  that  section  in  the 
preceding  problem. 

If  [{'units'}]  is  omitted  in  the  header  for  a  'New' 
section,  default  is  to  inches  and  pounds. 

j_.  If  data  are  added  to  a  preceding  section,  the  limitations 
on  number  of  lines  in  the  section, old  plus  added,  must  not 
exceed  those  specified  in  the  input  description  below. 

A  data  file  must  begin  with  a  Heading.  In  a  succession  of 
problems,  the  Heading  for  the  preceding  problem  may  be  re¬ 
placed  by  entering  the  new  heading  at  the  beginning  of  the 
new  problem  description.  Addition  to  a  previous  heading  is 
not  permitted. 

J_.  Comment  lines  may  be  inserted  in  the  input  file  by  enclos¬ 
ing  the  line,  following  the  line  number,  in  parentheses. 
Comment  lines  are  ignored,  e.g. , 

123l*|5  (THIS  LINE  IS  IGNORED). 
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Input  Description 


12.  Heading--One  (I)  to  four  (4)  lines  for  identifying  the  problem 
Line  contents 
LN  { 'heading'} 
b_.  Definition 

'heading'  =  any  alphanumeric  information  up  to  70  charac¬ 
ters  including  LN  and  any  imbedded  blanks. 

£.  Restriction:  If  a  'heading'  line  following  LN  begins  with 
a  combination  of  letters  which  are  permissible  abbrevia¬ 
tions  of  any  of  the  keywords  described  below,  the  'heading* 
must  begin  with  a  single  quote. 

13-  Beam  Cross  Section  Properties--Two  (2)  to  twenty-two  (22)  lines 
£.  Header--0ne  (l)  line 
Ml)  Contents 

LN  'Beam'  [{'New'  or  'Add'})  [{'units'}) 

£(2)  Definition 

'Beam'  =  section  title; 

[{'units'})  =  'j_nches'  or  'Feet'  and  'Pounds'  or 
'Kjps'  for  'New';  default  to  inches 
and  pounds  if  omitted  for  'New'. 
b_-  Data  Lines — One  (1)  to  twenty-one  (21)  lines 
Ml)  Contents 

LN  XI  X2  E  A I  SII  [  A2  SI2] 
b_(  2 )  Definitions 

XI  =  x-coordinate  at  beginning  of  distribution; 

X2  =  x-coordinate  at  end  of  distribution; 

E  =  modulus  of  elasticity  in  region  XI  to  X2; 

A)  =  cross  section  area  at  XI; 

SII  *  cross  section  moment  of  inertia  at  XI; 

A2.SI2  *  area  and  moment  of  inertia  at  X2. 

M3)  Oiscussion 

(a)  If  A2  and  SI2  are  omitted  area  and  moment  of 
inertia  are  assumed  to  be  constant  from  XI  to 
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(b)  If  A2  and  SI2  are  provided,  area  and  moment  of 
inertia  are  assumed  to  vary  linearly  from  Al 
and  SI1  at  XI  to  A2  (and  SI2  at  X2). 

(c)  When  several  distributions  are  described,  the 
minimum  of  all  XI  values  defines  the  left  end 
of  the  beam  and  the  maximum  of  al 1  X2  values 
defines  the  right  end  of  the  beam.  There  must 
be  positive,  nonzero  values  of  modulus  of  elas¬ 
ticity,  area,  and  moment  of  inertia  for  every 
x-coordinate  between  the  ends  of  the  beam. 

(d)  If  distributions  are  specified  for  overlapping 
regions,  values  within  the  overlap  are  cumula¬ 
tive. 

14.  Node  Spacing--One  (1)  to  twenty-two  (22)  lines 
£.  Header — One  (1)  line 
£(1)  Contents 

LN  1  NODe 1  [{'New'  or  'Md'}]  [{ 1  units1}] 

a^(2 )  Definitions 

1 NODe 1  =  section  title; 

[{'units'}]  =  'Jjiches'  or  'Feet'  for  'New' ;  de¬ 
fault  to  inches  if  omitted  for  'New1. 
b_.  Data  Line — One  (1)  to  twenty-one  (21)  lines 
b^(l)  Contents 

LN  XI  X2  HMAX 
b^(2)  Definitions 

XI, X2  =  x-coord i nates  at  beginning  and  end  of 
interval,  respectively; 

HMAX  =  maximum  distance  between  nodes  in  interval 
XI  to  X2 . 

£.  Discussion 

£(l)  Intervals  XI  to  X2  must  proceed  sequentially  from 
left  to  right.  Overlapping  intervals  or  gaps  be¬ 
tween  intervals  are  not  permitted. 

£(2)  XI  on  the  first  line  input  must  be  less  than  or 
equal  to  the  x-coordinate  at  the  left  end  of  the 
beam. 
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£(3)  X2  on  the  last  line  input  must  be  greater  than  or 
equal  to  the  x-coordinate  at  the  right  end  of  the 
beam. 

£(4)  Data  specified  beyond  the  ends  of  the  beam  are  ig¬ 
nored  . 

£(5)  Node  spacing  data  must  result  in  at  least  three  (3) 
nodes  (2  elements)  and  not  more  than  two-hundred-one 
(201)  nodes  (200  elements)  in  the  model. 

15.  Applied  Loads  —  Zero  (0)  or  two  (2)  to  twenty-two  (22)  lines. 
Entire  section  may  be  omitted 
a_.  Headet — One  (1)  line 
£(1)  Contents 

LN  1 LOads 1  [{'New1  or  'Add'}]  ({'units'}] 

£(2)  Definitions 

' LOads'  =  section  title; 

[{'units'}]  =  'Jnches'  or  'Feet'  and  'Pounds'  or 
'J<ips'  for  'New';  default  to  inches 
and  pounds  if  omitted. 

b^.  Data  Lines  for  Concentrated  Loads--0ne  (1)  line  for  each 
concentrated  load 
b^(l)  Contents 

LN  Concentrated'  XI  FX  FY  C 
b^2)  Definitions  for  concentrated  loads 
Concentrated'  =  keyword; 

XI  =  x-coordinate  at  point  of  concen¬ 
trated  load(s) : 

FX.FY  =  magnitudes  of  concentrated  x  and 
y  loads,  respectively; 

C  =  magnitude  of  concentrated  couple. 
£.  Data  Lines  for  Distributed  Loads--0ne  (1)  line  for  each 
d i st  r i but  ion 
£(1)  Contents 

LN  'Distributed'  'direction'  XI  Q1  X2  [Q2] 

£(  2 )  De  f  i  n  i  t  i  on  s 

'Distributed  =  keyword; 
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'  d i rection* 


*  'X'  for  axial  load,  or  1 Y *  for 
lateral  load; 

XI  =  x-coordinate  at  start  of  distribu¬ 
tion; 

X2  =  x-coordinate  at  end  of  distribution; 
Ql  =  magnitude  of  distributed  load  at  XI; 
[Q2]  =  magnitude  of  distributed  load  at  X2; 
assumed  to  be  equal  to  Ql  if  omitted 
( i  .e. ,  uni  form  load) . 

d.  Discussion 

d_(l)  Multiple  loads  specified  at  a  single  x-coordinate, 

either  concentrated  or  overlapping  distributions  are 
cumulat i ve. 

d_(2 )  Data  specified  beyond  the  ends  of  the  beam  are  ignor¬ 
ed. 

16.  Fixed  Supports--Zero  (0)  or  two  (2)  to  eleven  (II)  lines.  En¬ 
tire  section  may  be  omitted 
a.  Header- -One  (I)  line 
aO )  Contents 

LN  1 F I Xed 1  [{'New1  or  'Add'}]  [{'units'}] 

aj2)  Definitions 

1 F I Xed 1  =  section  title; 

[{'units'}]  =  'inches'  or  'Feet'  for  'New';  default 
to  inches  if  omitted. 
b_.  Data  Lines--0ne  (1)  to  ten  (10)  lines 
b^(l)  Contents 

LN  XI  { XD>  {YD}  { R} 
b^(2)  Definitions 

XI  ■  x-coordinate  at  point  of  support; 

{XD}  =  specified  X-displacement ;  enter  'Free'  if 
unspeci f ied; 

{YD}  *  specified  Y-displacement ;  enter  'Free'  if 
unspeci f ied; 

{R}  ■  specified  rotation  (in  radians);  enter 
'Free'  if  unspecified. 
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c. 


Discuss  ion 

£(1)  Multiple  values  of  a  single  displacement  component 
not  permitted. 

£(2)  Data  specified  beyond  the  ends  of  the  beam  are  ig- 
no  red . 

17.  Linear  Spring  Supports--Zero  (0)  or  two  (2)  to  twenty-two  (22) 
lines.  Entire  section  may  be  omitted 
£.  Headei — One  (1)  line 
a_(  1 )  Contents 

LN  1 L I Nea r 1  [{'New1  or  'Add'}]  [{'units'}] 

a(2)  Definitions 

' L  I  Near'  =  section  title; 

[{'units'}]  =  '£nches'  or  'Feet'  and  'Pounds'  or 
1  K_i  ps  ’  for  'New';  default  to  inches 
and  pounds  if  omitted. 

£.  Data  Lines  for  Concentrated  Linear  Springs--0ne  (1)  line 
for  each  concentrated  linear  spring 
£(1)  Contents 

LN  '£oncentrated '  XI  ANGLE  ST  SR 
£(2)  Definitions 

'£oncent rated '  =  keyword; 

XI  =  x-coordinate  at  point  of  attach¬ 
ment  ; 

ANGLE  =  angle  (in  degrees)  between  x-axis 
and  line  of  action  of  translation 
spring;  measured  positive  counter¬ 
clockwise  from  x-axis; 

ST  =  stiffness  of  translation  spring  in 
force  per  unit  displacement; 

SR  =  stiffness  of  rotation  spring  in 
moment  per  radian. 

£.  Data  Lines  for  Distributed  Linear  Springs--0ne  (I)  line  for 
each  distribution 
£(1)  Contents 

LN  'Distributed'  'direction'  XI  X2  A  B  C 


O 

A9 


£(2)  Definitions 

'Distributed'  =  keyword; 

'direction'  =  'X'  for  axial  spring  or  ' Y '  for 
lateral  spring; 

XI  =  x-coordinate  at  start  of  distribu¬ 
tion; 

X2  =  x-coordinate  at  end  of  distribution; 
A,B,C  =  spring  stiffness  coefficients  (see 
discussion  below). 

d_.  Discussion 

d^(l)  If  several  concentrated  linear  springs  are  attached 
to  a  single  point,  effects  are  cumulative, 
d (2 )  Stiffness  of  a  distributed  spring  is  assumed  to  vary 

C 

according  to  S  =  A  +  BZ  ,  where  Z  is  the  distance 
from  XI.  Units  of  A  and  B  are  force  per  unit  length 
per  unit  displacement;  C  is  a  dimensionless  constant. 
d_( 3 )  Overlapping  distributions  result  in  cumulative  effects 
in  the  overlap. 

d_(4)  Data  specified  beyond  the  end  of  the  beam  are  ignored. 
18.  Nonlinear  Spring  Supports --Zero  (0)  or  four  (4)  to  sixty-four 
(64)  lines.  Entire  section  may  be  omitted 
a.  Header--0ne  (1)  line 
£(1)  Contents 

LN  ' N0N1 inear'  [{'New'  or  'Add'}]  [{'units}] 
a (2)  Definitions 

' NON 1  inear'  =  section  title; 

[{'units'}]  =  'inches'  or  'Feet'  and  'Pounds'  or 

'.Kips';  default  to  inches  and  pounds 
i f  omi tted. 

b_.  Data  Lines  for  Concentrated  Nonlinear  Spr ings--Three  (3) 
lines  for  each  concentrated  spring 
b  ( 1 )  Line  1  Contents: 

LN  Concentrated'  XI  ANGLE  NPTS  DMUL  FMUL 
Line  2  Contents: 

LN  DEF(J)  DEF (2)  .  .  .  DEF(NPTS) 


Line  3  Contents: 

LN  FORCE ( 1 )  FORCE (2)  .  .  .  FORCE(NPTS) 
b_(2)  Definitions 

'£oncentrated 1  =  keyword 

XI  =  x-coordinate  at  point  of  attach¬ 
ment  ; 

ANGLE  =  angle  (in  degrees)  between  x-axis 
and  spring  line  of  action;  measur¬ 
ed  positive  counterclockwise  from 
positive  x-direction; 

NPTS  =  number  of  points  on  resisting 

force-deformation  curve;  minimum 
of  two  (2)  points  required;  maxi¬ 
mum  of  eight  (8)  points  permitted; 

DMUL  =  scale  factor  for  curve  deformation 
coordinate  values;  must  be  posi¬ 
tive  and  nonzero; 

FMUL  =  scale  factor  for  curve  resisting 
force  coordinate  values; 

DEF(  )  *  curve  deformation  coordinate  value; 

NPTS  values  on  one  line;  must  pro¬ 
ceed  sequentially  from  most  nega¬ 
tive  to  most  positive; 

F0RCE(  )  =  curve  resisting  force  coordinate 
value  corresponding  to  0 1 SP (  ). 

£.  Data  Lines  for  Distributed  Nonlinear  Springs — Minimum  of 
Two  (2)  groups  of  three  (3)  lines  for  each  distribution 
£(l)(a)  Group  1,  Line  1  Contents: 

LN  'Distributed1  'direction'  XI  NPTS  DMUL  FMUL 
Group  l.  Line  1  Contents: 

LN  DEF(l)  DEF (2)  .  .  .  DEF(NPTS) 

Group  I,  Line  3  Contents: 

LN  FORCE(I)  FORCE (2)  .  .  .  FORCE (NPTS) 
c_(  1 )  (b)  Definitions  for  Group  1: 

'Distributed'  *  keyword  indicating  beginning  of 
distribution; 


iz#**X,*  **JfW**^  "iWrtMiMfri'ltfi IK fi tk r ?i 1 1 


' d i rect ion1 


c  (2)  (a) 


c (2) (b) 


=  'X'  or  '£'  indicating  distribut¬ 
ed  spring  resists  axial  (x)  or 
lateral  (y)  displacements; 

XI  =  x-coordinate  at  beginning  of 
distribution; 

NPTS  =  number  of  points  provided  on  non¬ 
linear  resisting  force-displace¬ 
ment  curve  at  XI;  minimum  of  two 
(2)  points  required;  maximum  of 
eight  (8)  points  permitted; 

DMUL  =  scale  factor  for  curve  deforma¬ 
tion  coordinate  values;  positive, 
nonzero ; 

FMUL  =  scale  factor  for  curve  resisting 
force  coordinate  values; 

DEF(  )  =  deformation  coordinate  value; 


NPTS  values  on  one  line;  must 

proceed 

from  most  negative  to 

most  positive; 

FORCE (  ) 

=  resisting  distributed  force  cor¬ 

responding  to  D 1 SP (  ). 

Group 

2 ,  Line  1 

Contents : 

LN 

'control  ' 

X2  NPTS 

DMUL  FMUL 

Group 

2,  Line  2 

Contents  : 

LN 

DEF ( 1 )  DEF (2 )  .  .  . 

def(npts) 

Group 

2,  Line  3 

Contents : 

LN 

FORCE  (1 ) 

FORCE (2 )  . 

.  .  FORCE (NPTS) 

Def  in 

it  ions  for 

Group  2: 

'control'  =  'Continue'  or  '£nd';  '£ontinue'  indi¬ 
cates  distribution  continues  and  an¬ 
other  set  of  Group  2  lines  follows; 
'£nd'  indicates  this  is  last  curve 
in  this  distribution; 

X2  =  x-coordinate  for  this  nonlinear 
curve ; 

NPTS,  OMUL,  FMUL,  DEF (  ),  FORCE (  )  as  in  Group  1. 


At  2 


d_.  Discussion 

d^(l)  Final  deformation  and  resisting  force  curve  coordi¬ 
nates  are  products  DMUL’DEF(  )  and  FMUL*FORCE(  ), 
respectively. 

d^(2)  Resisting  force-deformation  curve  must  be  single 
valued  at  every  point. 

£(3)  Resisting  force-deformation  curve  is  assumed  to  vary 
linearly  between  input  points. 

d_(b)  Resisting  force  is  assumed  to  be  constant  at  first 

or  last  value  provided  on  curve  for  deformations  be¬ 
yond  the  range  of  deformation  coordinates  provided. 

d_(5)  Sever  concentrated  nonlinear  springs  may  be  attach¬ 
ed  a  single  point. 

d_( 6 )  Characteristics  of  distributed  nonlinear  springs  at 

intermediate  points  on  a  distribution  are  obtained  by 
linear  interpolation  between  adjacent  curves. 

dK7)  Overlapping  distributions  result  in  cumulative 
effects  in  the  overlap. 

c[(8)  Data  specified  beyond  the  ends  of  the  beam  are  ignor¬ 
ed. 

19.  Termi nat ion--0ne  (1)  line 
£.  Contes : - 

LN  1  F_l_Ni  sh 1  [{'Rerun'}] 

b_.  Definitions 

1 F I Ni sh 1  =  keyword  indicating  end  of  input  data  for  a 
problem.  Initiates  data  checking  and  solu¬ 
tion; 

[{'Rerun'}]  =  keyword;  if  omitted,  program  assumes  this 
is  end  of  data  file;  if  included,  after 
solution  of  preceding  problem,  program  will 
immediately  read  succeeding  data  as  input 
for  a  new  problem. 


Al  3 


'-itC'i 3.- 


Abbreviated  Input  Guide 


)ata  items  enclosed  in  brackets  [  ]  may  be  omitted.  Braces  {  }  en- 
ig  data  lists  indicate  choose  one;  arrow  indicates  default  if  item 
tted.) 

).  Heading — One  (1)  to  four  (4)  lines 
LN  'heading' 

[ LN  'heading'] 

[ LN  'heading'] 
l LN  'heading'] 

I.  Beam  Cross  Section  Propert i es--Two  (2)  to  twenty-two  (22)  lines 
Headei — One  (1)  line 

.  ,  r,  -*■  '  New' ,  ,  rf+'  Inches  1  Pounds  \  , 

LM  'Beam'  [{  K  ,^eet.  H  -Kips'  }1 

b.-  Data  Lines — One  (l)  to  twenty-one  (21)  lines 
LN  XI  X2  E  Al  SI1  [A2  SI2] 

!.  Node  Spacing--Two  (2)  to  twenty-two  (22)  lines 
a.  Header--One  (1)  line 

ln  -Nooe'  ir:g;:>]  ir:g2?s,n 

C-  Data  Lines--0ne  (1)  to  twenty-one  (21)  lines 
LN  XI  X2  HMAX 

1.  Applied  ..oads-'Zero  (0)  or  two  (2)  to  twenty-two  (22)  lines 
Header--0ne  (1)  line 

LN  'Loads'  [{  ,^d,}]  [{  ,|-eet,  H  -Kips-  11 

Data  Lines  for  Concentrated  Loads — One  (1)  line  for  each 
concentrated  load 

LN  'Concentrated'  XI  FX  FY  C 
£.  Data  Lines  for  Distributed  Loads--0ne  (1)  line  for  each 
distribution 


LN 


'Distributed'  {,*,}  XI  Q]  X2  [Q2] 


Fixed  Supports--Zero  (0)  or  two  (2)  to  eleven  (II)  lines 
a.  Header — One  (1)  line 


LN  'FlXed' 


r{-*-'New' } ,  [{^'!r 

I  A  A  A  I  ;  J  l  1  1C- 


'Add 


'  Jnches 1 
Feet ' 


)] 


AI4 


£•  Data  Lines--One  (1)  to  ten  (10)  lines 

LN  XI  {??  ,}  (™  ,}  {*  ,} 

’Free'  'Free'  'Free' 

25.  Linear  Spring  Supports--Zero  (0)  or  two  (2)  to  twenty-two  (22) 

1  ines 

£.  Header — One  (1)  line 

LN  UNear  [{  ,^d  }]  [{  ,peet,  H  ,^psI  H 

b  Data  Lines  for  Concentrated  Linear  Springs — One  (1)  line 
for  each  spring 

LN  'Concentrated'  XI  ANGLE  ST  SR 
£.  Data  Lines  for  Distributed  Linear  Springs--0ne  (1)  line  for 
each  distribution 

LN  'Distributed'  {,*,}  XI  X2  A  B  C 

26.  Nonlinear  Spring  Supports-- Zero  (0),  four  (4),  or  seven  (7)  to 
(64 )  I i nes 

a.  Header — One  (l)  line 


LN  'NONlinear'  [C  ,%£,)]  [C.g™*  }C 5}] 

£•  Data  Lines  for  Concentrated  Nonlinear  Springs — Three  (3) 
lines  for  each  spring 
Line  1  : 

LN  'Cor.  cent  rated1  XI  ANGLE  NPTS  DMUL  FMUL 


Line  2 : 


LN  DEF(l)  DEF (2)  .  .  .  DEF(NPTS) 

Line  3- 

LN  FORCE ( 1 )  FORCE (2)  .  .  .  FORCE(NPTS) 

c.  Data  Lines  for  Distributed  Nonlinear  Springs — Minimum  of 
two  (2)  groups  of  three  (3)  lines  for  each  distribution 
Group  1,  Line  1: 

LN  '£istributed'  {,*,}  XI  NPTS  DMUL  FMUL 
Group  1 ,  Line  2 : 

LN  DEF ( 1 )  DEF (2)  .  .  .  DEF(NPTS) 

Group  l ,  Line  3- 

LN  FORCE ( 1 )  FORCE (2)  .  .  .  FORCE(NPTS) 


Group  2,  Line  1 : 

LN  {!£°||?inUe'>  X2  NPTS  DMUL  FMUL 
Group  2,  Line  2: 

LN  DEF ( 1 )  DEF (2)  .  .  .  DEF(NPTS) 

Group  2,  Line  3: 

LN  FORCE(l)  FORCE (2)  .  .  .  FORCE(NPTS) 
27.  Termination 

LN  'FINish'  ['Rerun'] 


AI6 
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APPENDIX  B 
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EXAMPLE  SOLUTIONS 


Bl 


*****  input  file  for  example  1  ***** 


1000 
1010 
1020 
I  02.0 
1040 
3  050 
i.  060 
1070 
1.080 
.1  090 
1100 
1  1 1  0 
1  1.20 
1130 
1140 
1  1.50 
1160 
1170 
1  180 
1190 
1200 
1 2 1 0 
1720 
1  230 
1  2  40 
1250 
1260 
1270 
1280 
1 2  9  o 
1  300 
1310 
1370 
15  30 
1340 


EXAMPLE  1  --  FIXED  END  BEAM  ANALYSIS 
EFFECTS  OF  AXIAL  LOAD  ON  BENDING  EXCLUDED 
BEAM 

0  120  30. E6  4  40 
NODES  F 
-10  20  5 
LOADS  F  K 
D  Y  0  1.2  10 
r,  o  40  o  o 
FIXED  F 
0  FREE  0  0 
10  <;■  0  0 
FINISH  RERUN 

EXAMPLE  1A  --  SAME  AS  EXAMPLE  1 
EXCEPT  INCLUDE  EFFECTS  OF  AXIAL 
AXIAL  LOAD  =  40  NIPS 
FINISH  RERUN 

EXAMPLE  IB  —  INCREASE  AXIAL  LO 


BENDING 


EXAMPLE  IB  --- 
LOADS  ADD 
C  0  1960  0  0 
FINISH  RERUN 
EXAMPLE  1C  -- 
LOADS  ADD 
C  0  1000  0  0 
FINISH  RERUN 
EXAMPLE  ID  — 


AXIAL  LOAD  TO 


INCREASE  AXIAL  LOAD 


3000  KIPS 


INCREASE  AXIAL  LOAD 
ID  IS  SLIGHTLY  LESS 


TO  3289  KIPS 
THAN  BUCKLING  LOAD 


THIS  AXIAL  LOAD  IS  SLIGHTLY  LESS  THAN  BUCKLING  LOAD 
LOADS  ADD 
C  0  289  0  0 
FINISH  RERUN 

EXAMPLE  JE  --  INCREASE  AXIAL  LOAD  TO  3290  KIPS 

THIS  AXIAL  LOAD  IS  SLIGHTLY  GREATER  THAN  BUCKLING  LOAD 

LOADS  ADD 

C  0  1  0  0 

FINISH 


.  '  .  -* 

.  i  '  .  \ 


PROGRAM  CBEAMC  -  ANALYSIS  OF  BEAM-COLUMNS  WITH  NONLINEAR  SUPPORTS 
DATE.*  06/23/80  TIME*.  19556JOI 


ARE  INPUT  DATA  TO  BE  READ  FROM  TERMINAL  OR  FILE? 
ENTER  'TERMINAL'  OR  'FILE' 

I>F 

ENTER  INPUT  FILE  NAME  (6  CHARACTERS  MAXIMUM) 
I>EXAMP1 

INPUT  COMPLETE 

DO  YOU  UANT  INPUT  DATA  ECHOPRINTED  TO  YOUR 
TERMINAL*  TO  A  FILE*  TO  BOTH*  OR  NEITHER? 

ENTER  'TERMINAL'.  'FILE'.  'BOTH'.  OR  'NEITHER' 

I>T 


B4 


>  .  .  J-v.  \ 

■IS.;*  ;  •'  •  • 

■  .  .  .w- 

•  >—*»  ,  •. 


PROGRAM  CBEAMC  -  ANALYSIS  OF  BEAM-COLUMNS  WITH  NONLINEAR  SUPPORTS 
DATE!  06/23/00  TIMES  19S56S43 

I.  — INPUT  DATA 


I.— HEADING 

EXAMPLE  1  —  FIXED  END  BEAM  ANALYSIS 
EFFECTS  OF  AXIAL  LOAD  ON  BENDING  EXCLUDED 


2. --BEAM  CROSS  SECTION  DATA 


X-COORDINATE 
START  STOP 

< IN)  (IN) 

0.00  120.00 


MODULUS  OF 
ELASTICITY 
(PSD 
3.00E+07 


< - SECTION  PROPERTIES 


< - START - > 

AREA  INERTIA 

(SQIN)  (IN**4> 

4.00  40.00 


< - STOp - > 

AREA  INERTIA 

(SQIN)  (IN**4> 

4.00  40.00 


3.— NODE  SPACING  DATA 
X-COORDINATE 
START  STOP 

(FT)  (FT) 

-10.00  20.00 


MAXIMUM  NODE 
SPACING 
(FT) 

5.00 


4.—  LOAD  DATA 

4. A.— CONCENTRATED  LOADS 

CONCENTRATED  LOADS 


X-COORD 

X-LOAD 

Y-LOAD 

COUPLE 

(FT) 

(K) 

(K> 

(K-FT) 

0.00 

40.00 

0.00 

0.00 

4.B.— DISTRIBUTED 
LOAD 

LOADS 

DIRECT 

X-COORD 

LOAD 

X-COORD 

LOAD 

(FT) 

(K/FT) 

(FT) 

(K/FT) 

Y 

0.00 

1.20 

10.00 

1.20 

—FIXED  SUPPORT  DATA 
SUPPORT 

SPECIFIED  DISPLACEMENTS 

X-COORD 

X-DISP. 

Y-DISP. 

ROTATION 

(FT) 

(FT) 

(FT) 

(RAD) 

0.00 

FREE 

0.00 

0.00 

10.00 

0.00 

0.00 

0.00 

6.— LINEAR  SPRING  DATA 
NONE 


7,—  NONLINEAR  SPRING  DATA 
NONE 


DO  YOU  UANT  TO  EDIT  INPUT  DATA?  ENTER  'YES'  OR  'NO 


I>N 

DO  YOU  UANT  TO  CONTINUE  SOLUTION?  ENTER  'YES'  OR  'NO' 

I>Y 

DO  YOU  UANT  AXIAL  FORCE  EFFECTS  ON  SENDING  STIFFNESS  INCLUDED 
IN  THE  SOLUTION?  ENTER  'YES'  OR  'NO' 

I>N 


SOLUTION  COMPLETE 

DO  YOU  UANT  RESULTS  URITTEN  TO  YOUR  TERMINAL . 

TO  A  FILEt  OR  DOTH? 

ENTER  'TERMINAL'*  'FILE'.  OR  'DOTH' 

I>T 

ENTER  DESIRED  OUTPUT  UNITS 

'INCHES'  OR  'FEET'.  AND  'POUNDS'  OR  'KIPS'.  OR  'DEFAULT' 

I>D 


PROGRAM  CBEAMC  -  ANALYSIS  OF  BEAM-COLUMNS  UITH  NONLINEAR  SUPPORTS 
date:  06/23/eo  time:  i?:57:i4 


II. — SUMMARY  OF  RESULTS 
II. A. --HEADING 

EXAMPLE  1  —  FIXED  END  BEAM  ANALYSIS 
EFFECTS  OF  AXIAL  LOAD  ON  BENDING  EXCLUDED 


1 1.  B.— MAXIMA 


MAXIMUM 

X-COORD 

MAXIMUM 

X-COORD 

POSITIVE 

(IN) 

NEGATIVE 

(IN) 

AXIAL  DISPLACEMENT  <IN>  : 

4.000E-02 

0.00 

0. 

0.00 

LATERAL  DISPLACEMENT  (IN): 

4 . 500E-02 

AO.  00 

0. 

0.00 

ROTATION  (RAD)  1 

0. 

0.00 

0. 

0.00 

AXIAL  FORCE  <P>  : 

0. 

0.00 

-4.000E+04 

0.00 

SHEAR  (P)  : 

6.000E+03 

120.00 

-6.000E+03 

0.00 

BENDING  MOMENT  (P-IN)  : 

1 . 200E+03 

0.00 

-6 . 000E+04 

60.00 

II. C.— REACTIONS  AT  FIXED  SUPPORTS 

X-COORD  X-REACTION 

< IN)  <P> 

0.00  0. 

120.00  -4.00E+04 


Y-REACTION 

<P> 

-6.00E+03 

-6.00E+03 


NON-REACTION 

<P-IN> 

-1.20E+03 

1.20E+05 


ENTER  ONE  OF  FOLLOUING  TO  OBTAIN  COMPLETE  TABULATION  OF  RESULTS: 
'ALL'.  OR  RANGE  OF  X-COORDINATES  IN  INCHES  .  OR  'NONE'.  OR  'HELP' 

I>N 

DO  YOU  UANT  OUTPUT  UITH  DIFFERENT  UNITS?  ENTER  'YES'  OR  'NO' 

I>N 


OUTPUT  COMPLETE 
INPUT  COMPLETE 

DO  YOU  WANT  INPUT  DATA  ECHOPRINTED  TO  YOUR 
TERMINAL.  TO  A  FILE.  TO  BOTH»  OR  NEITHER? 

ENTER  ' TERMINAL ' »  'FILE'.  'BOTH'.  OR  'NEITHER' 

I>N 

00  YOU  NANT  TO  EDIT  INPUT  DATA?  ENTER  'YES'  OR  'NO' 

I>N 

DO  YOU  WANT  TO  CONTINUE  SOLUTION?  ENTER  'YES'  OR  'NO' 

I>Y 

DO  YOU  WANT  AXIAL  FORCE  EFFECTS  ON  BENDING  STIFFNESS  INCLUDED 
IN  THE  SOLUTION?  ENTER  'YES'  OR  'NO' 

I>Y 


SOLUTION  COMPLETE 

DO  YOU  WANT  RESULTS  WRITTEN  TO  YOUR  TERMINAL. 

TO  A  FILE.  OR  BOTH? 

ENTER  'TERMINAL'.  'FILE'.  OR  'BOTH' 

I>T 

ENTER  DESIRED  OUTPUT  UNITS 

'INCHES'  OR  'FEET'.  AND  'POUNDS'  OR  'KIPS'.  OR  'DEFAULT' 

I>D 


PROGRAM  CBEAMC  -  ANALYSIS  OF  BEAM-COLUMNS  WITH  NONLINEAR  SUPPORTS 
DATE!  06/23/80  TIME!  19:S8!01 


II.—  SUMMARY  OF  RESULTS 


II.  A. —HEADING 

EXAMPLE  1A  —  SAME  AS  EXAMPLE  1 

EXCEPT  INCLUDE  EFFECTS  OF  AXIAL  LOAD  ON  BENDING 

AXIAL  LOAD  «  40  KIPS 


1 1.  B.— MAXIMA 


AXIAL  DISPLACEMENT  <IN) 
LATERAL  DISPLACEMENT  (IN) 
ROTATION  (RAD) 

AXIAL  FORCE  (P) 

SHEAR  (P> 

BENDING  MOMENT  (P-IN) 


MAXIMUM 

X-COORD 

POSITIVE 

(IN) 

4 . 000E-02 

0.00 

4.333E-02 

60.00 

0. 

0.00 

0. 

0.00 

6 . 000E+03 

120.00 

1 . 210E+0S 

0.00 

MAXIMUM 

X-COORD 

NEGATIVE 

(IN) 

0. 

0.00 

0. 

0.00 

0. 

0.00 

■4 .  OOOE+04 

60.00 

6.000E+03 

0.00 

■6.0S5E+04 

60.00 

1 1.  C.—  REACTIONS  AT  FIXED  SUPPORTS 


X-COORD 

(IN) 

0.00 

120.00 


X— REACTION 
<P> 

0. 

-4.00E+04 


Y-REACTION 

(P> 

-6.00E+03 

-6.00E+03 


MOM-REACTION 

(P-IN) 

-1.21E+05 

1.21E+03 


ENTER  ONE  OF  FOLLOWING  TO  OBTAIN  COMPLETE  TABULATION  OF  RESULTS! 
'ALL'.  OR  RANGE  OF  X-COORDINATES  IN  INCHES  .  OR  'NONE'.  OR  'HELP' 

I>N 


I>u 


DO  YOU  WANT  OUTPUT  WITH  DIFFERENT  UNITS?  ENTER  'YES'  OR  'NO 


OUTPUT  COMPLETE 
INPUT  COMPLETE 

DO  YOU  WANT  INPUT  DATA  ECHOPRINTED  TO  YOUR 
TERMINAL.  TO  A  FILE.  TO  BOTH.  OR  NEITHER? 

ENTER  'TERMINAL'.  'FILE'.  'BOTH'.  OR  'NEITHER' 

I>N 

DO  YOU  WANT  TO  EDIT  INPUT  DATA?  ENTER  ' YE8'  OR  'NO' 

I>N 

DO  YOU  WANT  TO  CONTINUE  SOLUTION?  ENTER  'YES'  OR  'NO' 

I>Y 

DO  YOU  WANT  AXIAL  FORCE  EFFECTS  ON  BENDING  STIFFNESS  INCLUDED 
IN  THE  SOLUTION?  ENTER  'YES'  OR  'NO' 

I>Y 


SOLUTION  COMPLETE 

DO  YOU  WANT  RESULTS  WRITTEN  TO  YOUR  TERMINAL. 

TO  A  FILE.  OR  BOTH? 

ENTER  'TERMINAL'.  'FILE'.  OR  'BOTH' 

I>T 

ENTER  DESIRED  OUTPUT  UNITS 

'INCHES'  OR  'FEET'.  AND  'POUNDS'  OR  'KIPS'.  OR  'DEFAULT' 

I>D 


PROGRAM  CBEAMC  -  ANALYSIS  OF  BEAM-COLUMNS  WITH  NONLINEAR  SUPPORTS 
date:  oa/23/80  time:  i?:sa:«0 


1 1. —SUMMARY  OF  RESULTS 
II.  A. —HEADING 

EXAMPLE  IB  —  INCREASE  AXIAL  LOAD  TO  2000  KIPS 


II.  *.—  MAXIMA 


MAXIMUM 

X-COORD 

MAXIMUM 

X-COORD 

POSITIVE 

(IN) 

NEGATIVE 

(IN) 

AXIAL  DISPLACEMENT  (IN)  1 

2.000E+00 

0.00 

0. 

0.00 

LATERAL  DISPLACEMENT  (IN)1 

1.138E-01 

AO.  00 

0. 

0.00 

ROTATION  (RAD)  : 

0. 

0.00 

0. 

0.00 

AXIAL  FORCE  (P)  : 

0. 

0.00 

-2.000E+0A 

0.00 

SHEAR  (P>  : 

4.000E+03 

120.00 

-A.000E+03 

0.00 

BENDING  MOMENT  (P-IN)  ( 

2.373E+03 

0.00 

-1.703E+05 

AO.  00 

II. C 


—REACTIONS  AT  FIXED  SUPPORTS 


X-COORD 

(IN) 


X-REACTION 
<P> 


0.00 

120.00 


0. 

-2.00E+0A 


Y-RE ACTION 
<P> 

-A.00E+03 

-4.00E+03 


MOM-REACTION 
(P-IN) 
-2.37E+03 
2 • 37E+05 


ENTER  ONE  OF  FOLLOWING  TO  OBTAIN  COMPLETE  TABULATION  OF  RESULTS: 
'ALL'.  OR  RANGE  OF  X-COORDINATES  IN  INCHES  .  OR  'NONE'.  OR  'HELP' 

I>N 


DO  YOU  WANT  OUTPUT  WITH  DIFFERENT  UNIT8?  ENTER  'YES'  OR  'NO' 

I>N 


B8 


OUTPUT  COMPLETE 
INPUT  COMPLETE 

DO  YOU  WANT  INPUT  DATA  ECHOPRINTED  TO  YOUR 
TERMINAL*  TO  A  PILE*  TO  BOTH.  OR  NEITHER? 

ENTER  'TERMINAL' >  'PILE'.  'BOTH'.  OR  'NEITHER' 

I>N 

DO  YOU  WANT  TO  EDIT  INPUT  DATA?  ENTER  ' YES'  OR  'NO' 

I>N 

DO  YOU  WANT  TO  CONTINUE  SOLUTION?  ENTER  'YES'  OR  'NO' 

I>Y 

DO  YOU  WANT  AXIAL  FORCE  EFFECTS  ON  BENDING  STIFFNESS  INCLUDED 
IN  THE  SOLUTION?  ENTER  'YES'  OR  'NO' 

I>V 


SOLUTION  COMPLETE 

DO  YOU  WANT  RESULTS  WRITTEN  TO  YOUR  TERMINAL. 

TO  A  FILE.  OR  BOTH? 

ENTER  'TERMINAL'.  'FILE'.  OR  'BOTH' 

I>T 

ENTER  DESIRED  OUTPUT  UNITS 

'INCHES'  OR  'FEET'.  AND  'POUNDS'  OR  'KIPS'.  OR  'DEFAULT' 

I>D 


PROGRAM  CBEAMC  -  ANALYSIS  OF  BEAM-COLUMNS  WITH  NONLINEAR  SUPPORTS 
DATE}  06/23/80  TIME*  19559S34 


I I. --SUMMARY  OF  RESULTS 


1 1.  A.—  HEADING 

EXAMPLE  1C  —  INCREASE  AXIAL  LOAD  TO  3000  KIPS 


II.  B.  —MAXIMA 


AXIAL  DISPLACEMENT  (IN) 
LATERAL  DISPLACEMENT  (IN) 
ROTATION  (RAD) 

AXIAL  FORCE  (P> 

SHEAR  (P) 

BENDINO  MOMENT  (P-IN) 


MAXIMUM 

X-COORD 

POSITIVE 

(IN) 

3 . OOOE+OO 

0.00 

5.041E-01 

60.00 

0. 

0.00 

0. 

0.00 

6.000E+03 

120.00 

8.818E+03 

0.00 

MAXIMUM 

X-COORD 

NEGATIVE 

(IN) 

0. 

0.00 

0. 

0.00 

0. 

0.00 

3.000E+06 

0.00 

6.000E+03 

0.00 

8.103E+03 

60.00 

1 1. C. —REACT IONS  AT  FIXED  SUPPORTS 


T-COORD 

(IN) 

0.00 

120.00 


X-REACTION 

(P) 

0. 

-3.00E+06 


Y-REACTION 

(P> 

-6.00E+03 
— 6.00E+03 


MOM-REACTION 

(P-IN) 

-8.B2E+05 

B.82E+03 


ENTER  ONE  OF  FOLLOWING  TO  OBTAIN  COMPLETE  TABULATION  OF  RESULTS: 
'ALL'.  OR  RANGE  OF  X-COORDINATES  IN  INCHES  .  OR  'NONE'.  OR  'HELP' 

I>N 


DO  YOU  WANT  OUTPUT  WITH  DIFFERENT  UNITS?  ENTER  'YES'  OR  'NO' 

I>N 


B9 


OUTPUT  COMPLETE 
INPUT  COMPLETE 

DO  YOU  WANT  INPUT  DATA  ECHOPRINTED  TO  YOUR 
TERMINAL*  TO  A  PILE*  TO  BOTH*  OR  NEITHER? 

ENTER  ' TERMINAL '  *  'FILE'.  'BOTH'.  OR  'NEITHER' 

I>N 

DO  YOU  WANT  TO  EDIT  INPUT  DATA?  ENTER  'YES'  OR  'NO' 

I>N 

DO  YOU  WANT  TO  CONTINUE  SOLUTION?  ENTER  'YES'  OR  'NO' 

I>Y 

DO  YOU  WANT  AXIAL  FORCE  EFFECTS  ON  BENDING  STIFFNESS  INCLUDED 
IN  THE  SOLUTION?  ENTER  'YES'  OR  'NO' 

I>Y 


SOLUTION  COMPLETE 

DO  YOU  WANT  RESULTS  WRITTEN  TO  YOUR  TERMINAL* 

TO  A  FILE.  OR  BOTH? 

ENTER  'TERMINAL'.  'FILE'.  OR  'BOTH' 

I>T 

ENTER  DESIRED  OUTPUT  UNITS 

'INCHES'  OR  'FEET'.  AND  'POUNDS'  OR  'KIPS'.  OR  'DEFAULT' 

I>D 


PROGRAM  CBEAMC  -  ANALYSIS  OF  BEAM-COLUMNS  WITH  NONLINEAR  SUPPORTS 

date:  oa/23/bo  time:  20:00:22 


II.—  SUMMARY  OF  RESULTS 
II. A. —HEADING 

EXAMPLE  ID  —  INCREASE  AXIAL  LOAD  TO  3289  KIPS 
THIS  AXIAL  LOAD  IS  SLIGHTLY  LESS  THAN  BUCKLING  LOAD 


II  .B.— MAXIMA 


AXIAL  DISPLACEMENT  <IN>  : 
LATERAL  DISPLACEMENT  <IN>: 
ROTATION  (RAD)  : 
AXIAL  FORCE  (P)  : 
SHEAR  (P)  ; 
BENDING  MOMENT  (P-IN)  : 


MAXIMUM 

X-COORD 

POSITIVE 

(IN) 

3.289E+00 

0.00 

1.6BIE+02 

AO.  00 

0. 

0.00 

0. 

0.00 

A.000E+03 

120.00 

2.7A3E+08 

0.00 

MAXIMUM 

X-COORD 

NEGATIVE 

(IN) 

0. 

0.00 

0. 

0.00 

0. 

0.00 

3.289E+0A 

0.00 

•A . 00OE+03 

0.00 

2.7A4E+0B 

AO.  00 

1 1.  C.— REACTIONS  AT  FIXED  SUPPORTS 


X-COORD 

(IN) 

0.00 

120.00 


X-REACTION 

(P> 

0. 

-3.29E+0A 


Y-REACTION 

(P) 

-6.00E+03 

-A.00E+03 


MOM-REACTION 

(P-IN) 

-2.77E+08 

2.77E+08 


ENTER  ONE  OF  FOLLOWING  TO  OBTAIN  COMPLETE  TABULATION  OF  RESULTS. 
'ALL'.  OR  RANGE  OF  X-COORDINATES  IN  INCHES  ,  OR  'NONE'.  OR  'HELP' 

I>N 


DO  YOU  WANT  OUTPUT  WITH  DIFFERENT  UNITS?  ENTER  'YES'  OR  'NO' 

I>N 


BIO 


OUTPUT  COMPLETE 
INPUT  COMPLETE 

DO  YOU  WANT  INPUT  DATA  ECHOPRINTED  TO  YOUR 
TERMINAL,  TO  A  FILE »  TO  BOTH,  OR  NEITHER? 
ENTER  'TERMINAL',  'FILE',  'BOTH',  OR  'NEITHER' 
T 


HI  I 


1C  -  ANALYSIS  OF  BEAM-COLUMNS  WITH  NONLINEAR  SUPPORTS 
10  TIMES  20S00S47 

A 


•CREASE  AXIAL  LOAD  TO  3290  KIPS 

•AD  IS  SLIGHTLY  GREATER  THAN  BUCKLING  LOAD 


IS  SECTION  DATA 


< - - 

- SECTION 

PROPERTIES 

MODULUS  OF 

< - 

START - > 

< - 

STOP - > 

ELASTICITY 

AREA 

INERTIA 

AREA 

INERTIA 

(PSI) 

(SOIN) 

( IN. .4 ) 

( SGIN ) 

(IN**4> 

3.00E+07 

4.00 

40.00 

4.00 

40.00 

ING  DATA 

NATE  MAXIMUM  NODE 


STOP 

SPACING 

(FT) 

(FT) 

20.00 

5.00 

NTRATED 

LOADS 

CONCENTRATED 

LOADS 

OORD 

X-LOAD 

Y-LOAD 

COUPLE 

FT) 

(K) 

(K) 

(K-FT) 

.00 

40.00 

0.00 

0.00 

.00 

I960. 00 

0.00 

0.00 

.00 

1000.00 

0.00 

0.00 

.00 

289.00 

0.00 

0.00 

.00 

1.00 

0.00 

0.00 

IBUTED 

LOADS 

AD 

ECT 

X-COORD 

LOAD 

X-COORD 

LOAD 

(FT) 

(K/FT) 

(FT) 

(K/FT) 

1 

0.00 

1.20 

10.00 

1.20 

•ORT  DATA 


SPECIFIED  DISPLACEMENTS 
X-DISP.  Y-DISP.  ROTATION 

<  FT )  (FT)  (RAD) 

FREE  0.00  0.00 

0.00  0.00  0.00 


IING  DATA 


SPRING  DATA 


DO  YOU  WANT  TO  EDIT  INPUT  DATA?  ENTER  'YES'  OR  NO 


I  >N 

DO  YOU  UANT  TO  CONTINUE  SOLUTION?  ENTER  'YES'  OR  'NO' 

I>Y 

DO  YOU  UANT  AXIAL  FORCE  EFFECTS  ON  BENDING  STIFFNESS  INCLUDED 
IN  THE  SOLUTION’  ENTER  'YES'  OR  'NO' 

I>Y 


SOLUTION  COMPLETE 

DO  YOU  WANT  RESULTS  WRITTEN  TO  YOUR  TERMINAL. 

TO  A  FILE.  OR  BOTH? 

ENTER  'TERMINAL'.  'FILE'.  OR  'BOTH' 

I>T 

ENTER  DESIRED  OUTPUT  UNITS 

'INCHES'  OR  'FEET'.  AND  'POUNDS'  OR  'KIPS'.  OR  'DEFAULT' 


PROGRAM  CBEAMC  -  ANALYSIS  OF  BEAM-COLUMNS  WITH  NONLINEAR  SUPPORTS 
DATE!  06/23/80  TIME!  20:0i:i2 


II.— SUMMARY  OF  RESULTS 
II. A. —HEADING 

EXAMPLE  IE  INCREASE  AXIAL  LOAD  TO  3290  KIPS 

THIS  AXIAL  LOAD  IS  SLIGHTLY  GREATER  THAN  BUCKLING  LOAD 


II. B.— MAXIMA 


MAXIMUM 

X-COORD 

MAXIMUM 

X-COORD 

POSITIVE 

(IN) 

NEGATIVE 

(IN) 

AXIAL  DISPLACEMENT  (IN) 

3.290E+00 

0.00 

0. 

0.00 

LATERAL  DISPLACEMENT  (IN) 

0. 

0.00 

-1 . 106E+03 

60.00 

ROTATION  (RAD) 

0. 

0.00 

0. 

0.00 

AXIAL  FORCE  <P> 

0. 

0.00 

-3.290E+06 

0.00 

SHEAR  <P> 

6 . OOOE+03 

120.00 

-6. 000E+03 

0.00 

BENDING  MOMENT  (P-IN) 

1.820E+09 

60.00 

-1 .820E+09 

0.00 

II. C.—  REACTIONS  AT  FIXED 

SUPPORTS 

X-COORD 

X-REACTION 

Y-REACTION  MOM 

-REACTION 

(IN) 

<P> 

(P) 

(P-IN) 

0.00 

0. 

-6.00E+03 

1 .82E+09 

120.00 

-3.29E+06 

-6.00E+03 

-1 .82E+09 

ENTER  ONE  OF  FOLLOUING  TO  OBTAIN  COMPLETE  TABULATION  OF  RESULTS: 

'ALL'.  OR  RANGE  OF  X-COORDINATES  IN  INCHES  .  OR  'NONE'.  OR  'HELP' 

I>N 

DO  YOU  UANT  OUTPUT  WITH  DIFFERENT  UNITS?  ENTER  'YES'  OR  'NO' 

I>N 

OUTPUT  COMPLETE 

DO  YOU  UANT  TO  EDIT  INPUT  DATA  FOR  THE  PROBLEM 
JUST  COMPLETED?  ENTER  'YES'  OR  'NO' 

I>N 

DO  YOU  UANT  TO  MAKE  ANOTHER  RUN?  ENTER  'YES'  OR  'NO' 

I>N 

**********  NORMAL  TERMINATION  ********** 


HI  3 


'.dOMr.1:  -  — 


i  m 


*****  input  file  for  example  2  ***** 


1000  EXAMPLE  2  --  SIMPLE  BEAM  ON  UNIFORM  ELASTIC  FOUNDATION 
1010  ILLUSTRATE  EFFECT  OF  NODE  SPACING  ON  SOLUTION 
1020  'NODES  AT  5  FT. 

1030  BEAM 

1040  0  120  30. E6  4  40 
1050  NODES  F 
1060  0  10  5 
1070  LOADS  F  K 
1080  D  Y  0  1.2  10 
:  «9(>  FIXED  F 
1 100  O  FREE  0  FREE 
1.110  10  0  0  FREE 
1120  I  INEAR 

1130  D  Y  0  120  100  0  0 

I. 140  FINISH  RERUN 

1150  EXAMPLE  2A  --  SAME  AS  EXAMPLE  2  EXCEPT  NODES  AT  1.25  FT. 
1160  NODES  N  F 
1170  0  10  1.25 

I I. 80  FINISH 


131  ‘j 


PROGRAM  CBEAMC  -  ANALYSIS  OF  BEAM-COLUMNS  WITH  NONLINEAR  SUPPORTS 
DATE!  06/23/80  TIME:  20:22:18 

I.  — INPUT  DATA 

1.—  HEADING 

EXAMPLE  2  —  SIMPLE  BEAM  ON  UNIFORM  ELASTIC  FOUNDATION 
ILLUSTRATE  EFFECT  OF  NODE  SPACING  ON  SOLUTION 
'NODES  AT  5  FT. 


2.—  BEAM  CROSS  SECTION  DATA 


X-COORDINATE 

MODULUS  OF 

- SECTION 

START - > 

PROPERTIES 

< - 

STOP - > 

START 

STOP 

ELASTICITY 

AREA 

INERTIA 

AREA 

INERTIA 

(IN) 

(IN) 

(PS1) 

( SOIN) 

(IN**4) 

(SOIN) 

(IN**4> 

0.00 

120.00 

3.00E+07 

4.00 

40.00 

4.00 

40.00 

3.— NODE  SPACING  DATA 
X-COORDINATE 
START  STOP 

<FT>  (FT) 

0.00  10.00 


MAXIMUM  NODE 
SPACING 
(FT) 

S.00 


4.  —LOAD  DATA 

4. A.— CONCENTRATED  LOADS 
NONE 


4.  B.~ DISTRIBUTED  LOADS 


LOAD 

DIRECT 

< - STAF 

X-COORD 

!T - > 

LOAD 

< - STOP- 

X-COORD 

- > 

LOAD 

(FT) 

(K/FT) 

(FT) 

(K/FT) 

Y 

0.00 

1.20 

10.00 

1.20 

FIXED  SUPPORT 
SUPPORT 

DATA 

SPECIFIED  DISPLACEMENTS 

X-COORD 

X-DISP. 

Y-DISP . 

ROTATION 

(FT) 

(FT) 

(FT) 

(RAD) 

0.00 

FREE 

0.00 

FREE 

10.00 

0.00 

0.00 

FREE 

LINEAR  SPRING 

DATA 

6. A. --CONCENTRATED  LINEAR  SPRINGS 
NONE 


6.B. — DISTRIBUTED  LINEAR  SPRINGS 
SPRING  X -COORD 

DIRECTION  START  STOP 

(IN)  (IN) 

Y  0.00  120.00 


SPRING  STIFFNESS  COEFFICIENTS 
ABC 
(P/I/I)  (P/I/I) 

100.00  0.00  0.00 


7. --NONLINEAR  SPRING  DATA 
NONE 


\  5 


PROGRAM  CBEAMC 
DATE:  06/23/80 


ANALYSIS  OF  BEAN-COLUHNS  UITH  NONLINEAR  SUPPORTS 

tine:  20:22:44 


II. 


-SUMMARY  OF  RESULTS 


II.  A. —HEADING 

EXAMPLE  2  —  SIMPLE  BEAM  ON  UNIFORM  ELASTIC  FOUNDATION 
ILLUSTRATE  EFFECT  OF  NODE  SPACING  ON  SOLUTION 
'NODES  AT  5  FT. 


1 1.  B.— MAXIMA 


MAXIMUM 

X-COORD  MAXIMUM  X 

-COORD 

POSITIVE 

(IN)  NEGATIVE 

(IN) 

AXIAL  DISPLACEMENT  (IN) 

0. 

0.00  0. 

0.00 

LATERAL  DISPLACEMENT  (IN) 

1.912E-01 

60.00  0. 

0.00 

ROTATION  (RAD) 

S.116E-03 

0.00  -5.116E-03 

120.00 

AXIAL  FORCE  (P) 

0. 

0.00  0. 

0.00 

SHEAR  (P> 

5.273E+03 

120.00  -5.273E+03 

0.00 

BENDING  MOMENT  (P-IN) 

0. 

0.00  -1.522E+05 

60.00 

1 1. C.— REACTIONS  AT  FIXED 

SUPPORTS 

X-COORD 

X-REACTION  Y-REACTION  MOM-REACTION 

(IN) 

(P> 

(P>  (P- 

IN) 

0.00 

0. 

-5.27E+03  0. 

120.00 

0. 

-5.27E+03  0. 

PROGRAM  CBEAMC  -  ANALYSIS 

OF 

BEAM-COLUMNS 

UITH  NONLINEAR  SUPPORTS 

DATE:  06/23/80 

time:  20:22:44 

III.— COMPLETE  RESULTS 


III. A. —HEADING 

EXAMPLE  2  —  SIMPLE  BEAM  ON  UNIFORM  ELASTIC  FOUNDATION 
ILLUSTRATE  EFFECT  OF  NODE  SPACING  ON  SOLUTION 
'NODES  AT  3  FT. 


1 1 1.  B.—  DISPLACEMENTS  AND  INTERNAL  FORCES 


X-COORD 

(IN) 

0.00 

60.00 

120.00 


AXIAL 

(IN) 

0. 

0. 

0. 


-DISPLACEMENTS - > 

LATERAL  ROTATION 
(IN)  (RAD) 

0.  5.116E-03 

1.912E-01  0. 

0.  -5. 116E-03 


AXIAL 

(P) 


-INTERNAL  FORCES - 

SHEAR  MOMENT 

(P) 

-5.273E+03 

0. 

5 . 273E+03 


(P-IN) 


1 .522E+05 

0. 


III. C.~ FORCES  IN  DISTRIBUTED  LINEAR  SPRINGS 
DISTRIBUTED  SPRING  FORCES 


X-COORD 

(IN) 

0.00 

60.00 

120.00 


AXIAL 

(P/IN) 

0. 

0. 

0. 


LATERAL 

(P/IN) 

0. 

-1.912E+01 

0. 


B I  7 


1 


.  ..  *  V  *  . 


PROGRAM  CBEAMC  -  ANALYSIS  OF  BEAH-COLUHNS  WITH  NONLINEAR  SUPPORTS 
date:  06/23/80  tine:  20:23:28 


I I. --SUNNARY  OF  RESULTS 
I I. A. --HEADING 

EXAHPLE  2A  —  SANE  AS  EXAMPLE  2  EXCEPT  NODES  AT  1.25  FT. 


II. B.— MAXIMA 

MAXIMUM 

X-COORD 

MAXIMUM 

X-COORD 

POSITIVE 

(IN) 

NEGATIVE 

(IN) 

AXIAL  DISPLACEMENT  (IN)  :  0. 

0.00 

0. 

0.00 

LATERAL  DISPLACEMENT  (IN) 

1.910E-01 

60.00 

0. 

0.00 

ROTATION  (RAD) 

5. 109E-03 

0.00 

-5. 109E-03 

120.00 

AXIAL  FORCE  (P> 

0. 

0.00 

0. 

0.00 

SHEAR  (P) 

5.266E+03 

120.00 

-5. 266E+03 

0.00 

BENDING  MOMENT  (P-IN) 

0. 

0.00 

-1.520E+05 

60.00 

II.C. — REACTIONS  AT  FIXED 
X-COORD 

SUPPORTS 

X-REACTION 

Y-REACTION 

MOM-REACTION 

(IN) 

(P) 

(P> 

(P-IN) 

0.00 

0. 

-3.27E+03 

0. 

120.00 

0. 

-5.27E+03 

0. 

i 


; 


PROGRAM  CBEAMC  -  ANALYSIS  OF  BEAM-COLUMNS  WITH  NONLINEAR  SUPPORTS 
DATE!  06/23/80  TIME!  20:23 .*28 


I I I. --COMPLETE  RESULTS 
1 1 1.  A.— HEADING 

EXAMPLE  2A  —  SAME  AS  EXAMPLE  2  EXCEPT  NODES  AT  1.25  FT. 


I I I. B. --DISPLACEMENTS  AND  INTERNAL  FORCES 


:-C00RD 

< - 

AXIAL 

— DISPLACEHEN1 
LATERAL 

rs - > 

ROTATION 

< - 

AXIAL 

-INTERNAL  FORC 
SHEAR 

:es - > 

MOMENT 

(IN) 

(IN) 

(IN) 

(RAD) 

(P> 

(P> 

(P-IN) 

0.00 

0. 

0. 

5.109E-03 

0. 

-5.266E+03 

0. 

15.00 

0. 

7 . 434E-02, 

4.661E-03 

0. 

-3 . 823E+03 

-6.802E+04 

30.00 

0. 

1 . 362E-01 

3.495E-03 

0. 

-2 • 483E+03 

-t  * 152E  +  05 

45.00 

0. 

1 . 769E-01 

1 . 862E-03 

0. 

-1.221E+03 

-1.429E+05 

60.00 

0. 

1.910E-01 

0. 

0. 

0. 

-1.520E+05 

75.00 

0. 

1.769E-01 

-1.862E-03 

0. 

1 . 221E+03 

-1.429E+05 

VO.  00 

0. 

1 . 362E-01 

-3 . 495E-03 

0. 

2 • 483E+03 

-1.152E+05 

105.00 

0. 

7.434E-02 

-4.661E-03 

0. 

3 . 823E+03 

-6.802E+04 

120.00 

0. 

0. 

-5. 109E-03 

0. 

5.266E+03 

0. 

III. C.~ FORCES  IN  DISTRIBUTED  LINEAR  SPRINGS 


DISTRIBUTED 

SPRING  FORCES 

-COORD 

AXIAL 

LATERAL 

(IN) 

(P/IN) 

(P/IN) 

0.00 

0. 

0. 

15.00 

0. 

-7.434E+00 

30.00 

0. 

-1 . 362E+01 

45.00 

0. 

-1 • 769E401 

60.00 

0. 

-1 .910E+01 

75.00 

0. 

-1 • 769E+01 

90.00 

0. 

-1 • 362E+01 

105.00 

0. 

-7.434E+00 

120.00 

0. 

0. 

I 


O 


BI9 


♦swa 


Cross  Section:  E  =  4.3  x  10  ps i 
A  =  100  sq  in. 

I  =  833-33  in.4 

Figure  B3.  Example  3:  System  for  Pile  Head 
Stiffness  Matrix 


*****  input  file  for  example  3  ***** 


1000 
1010 
1020 
1030 
1040 
1050 
10 &0 
1070 

i  oao 

10  70 
1 100 
1110 
1  120 
11.30 
1140 
.1  130 
1160 
1170 

1 1  Mu 
!  !  90 
1200 
1210 
1  220 
1  2  30 
1  2  40 


EXAMPLE  3  ---  PILE  STIFFNESS  MATRIX 
AXIAL  FORCE  EFFECTS  OF  RENDING  EXCLUDED 
' NODES  AT  5FT. 

AXIAL  AND  LATERAL  TERMS 
BEAM 

0  1200  4.3EA  100  833.33 

NODES  F 

0  100  5 

FIXED  l 

Olio 

LINEAR 

D  Y  0  1200  0  10  2 
D  X  0  1200  0  10  2 
FINISH  RERUN 
EXAMPLE  3  --  CONTINUED 
AXIAL  AND  LATERAL  TERMS 
'NODES  AT  1  FT. 

NODES  NEW  F 
0  100  1 
FINISH  RERUN 
EXAMPLE  3  --  CONTINUED 
MOMENT  TERMS 
FIXED  NEW 
0  0  0  1 
FINISH 
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PROGRAM  CBEAHC  -  ANALYSIS  OF  BEAM-COLUMNS  WITH  NONLINEAR  SUPPORTS 
DATES  06/24/80  TINES  13S23S08 


I. --INPUT  DATA 
1.—  HEADING 

EXAMPLE  3  —  PILE  STIFFNESS  MATRIX 
AXIAL  FORCE  EFFECTS  ON  BENDING  EXCLUDED 
'NODES  AT  5  FT. 

AXIAL  AND  LATERAL  TERMS 


2. --BEAM  CROSS  SECTION  DATA 

<-COORDINATE  MODULUS  OF 

START  STOP  ELASTICITY 

<  IN)  (IN)  (PSD 

0.00  1200.00  4.30E+06 


<  - SECTION  PROPERTIES - > 

<  - START - >  < - STOP - > 

AREA  INERTIA  AREA  INERTIA 

(SQIN)  (IN**4>  <  SQIN)  (IN**4) 

100.00  833.33  100.00  833.33 


3.— NODE  SPACING  DATA 
X-COORDINATE 
START  STOP 

(FT)  (FT) 

0.00  100.00 


MAXIMUM  NODE 
SPACING 
(FT) 

5.00 


4. — LOAD  DATA 
NONE 


5.—  FIXED  SUPPORT  DATA 
SUPPORT 
X-COORD 
(IN) 

0.00 


SPECIFIED  DISPLACEMENTS 
X-DISP.  Y-DISP •  ROTATION 

(IN)  (IN)  (RAD) 

1.00  1.00  0.00 


6. --LINEAR  SPRING  DATA 

6. A.— CONCENTRATED  LINEAR  SPRINGS 
NONE 


6. B.— DISTRIBUTED  LINEAR  SPRINGS 
SPRING  X-COORD 

DIRECTION  START  STOP 

(IN)  (IN) 

Y  0.00  1200.^4 

X  0 . 00  1 200  ■  *0 


SPRING  STIFFNESS  COEFFICIENTS 
ABC 
(P/I/I)  (P/I/I) 

0.00  10.00  2.00 

0.00  10.00  2.00 


7. --NONLINEAR  SPRING  DATA 
NONE 
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PROGRAM  C8EAMC  -  ANALYSIS  OF  BEAM-COLUMNS  UITH  NONLINEAR  SUPPORTS 
DATE:  06/24/80  TIME!  I4:02!02 


II.—  SUMMARY  OF  RESULTS 
II. A. —HEADING 

EXAMPLE  3  --  PILE  STIFFNESS  MATRIX 
AXIAL  FORCE  EFFECTS  OF  BENDING  EXCLUDED 
'NODES  AT  5FT. 

AXIAL  AND  LATERAL  TERMS 


II .B. --MAXIMA 


AXIAL  DISPLACEMENT  (IN) 
LATERAL  DISPLACEMENT  (IN) 
ROTATION  (RAD) 

AXIAL  FORCE  (P) 

SHEAR  (P) 

BENDING  MOMENT  (P-IN) 


MAXIMUM  X-COORD 

POSITIVE  (IN) 

l.OOOE+OO  0.00 

l.OOOE+OO  0.00 

3 . 749E-05  120.00 

2.234E+02  300.00 

3.537E+05  0.00 

1.690E+06  60.00 


MAXIMUM  X-COORD 

NEGATIVE  (IN) 

-8 . 297E-06  300.00 

-1.547E-04  180.00 

-6.680E-03  60.00 

-3 . 892E+06  0.00 

-5 . 108E+04  60.00 

-7 . 331 E+06  0.00 


II. C. --REACTIONS  AT  FIXED  SUPPORTS 


X-COORD 

(IN) 

0.00 


X-REACT I  ON 
<P> 

3.89E+06 


Y -  RE ACT  ION 
(P) 

3.S4E+05 


MOM-REACTION 

(P-IN) 

7.33E+06 


PROGRAM  CBFAMC  -  ANALYSIS  OF  BEAH-COLUMNS  UITH  NONLINEAR  SUPPORTS 
date:  06/24/80  time:  14:02:02 


III.— COMPLETE  RESULTS 
III .A. --HEADING 

EXAMPLE  3  --  FILE  STIFFNESS  MATRIX 
AXIAL  FORCE  EFFECTS  OF  BENDING  EXCLUDED 
NODES  AT  5FT. 

AXIAL  AND  LATERAL  TERMS 


III. B. --DISPLACEMENTS  AND  INTERNAL  FORCES 


DISPLACEMENTS . . >  , . INTERNAL  FORCES 


X-COORD 

AXIAL 

LATERAL 

ROTATION 

AXIAL 

SHEAR 

MOMENT 

(  IN) 

<  IN) 

<  IN) 

(RAO) 

(P> 

(P) 

(P-IN) 

0.00 

l.OOOE+OO 

l.OOOE+OO 

0. 

-3 . 892E+06 

3.337E+03 

-7 . 3T 1E+V 

60.00 

4.944E-01 

4.9S9E-02 

-6 . 680E-03 

-3. 176E+06 

-5 . 108E  +  04 

1 .690E+C 

120.00 

1 .5B0E-01 

7.693E-04 

3.749E-05 

-1 .596E+06 

4 . 434E+03 

-S.210E+C 

III. C.— FORCES  IN  DISTRIBUTED  LINEAR  SPRINGS 
DISTRIBUTED  SPRING  FORCES 
X-COORD  AXIAL  LATERAL 

(IN)  (P/IN)  (P/IN) 

0.00  0.  0. 

60.00  -1.780E+04  -1.783E+03 

120.00  -2.276E+04  -1.108E+02 
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PROGRAM  CBEAMC 
DATE:  06/24/80 


-  ANALYSIS  OF  BEAN-COLUMNS  WITH  NONLINEAR  SUPPORTS 

time:  14:02:53 


II. --SUMMARY  OF  RESULTS 

II .A. --HEADING 
EXAMPLE  3  --  CONTINUED 
AXIAL  AND  LATERAL  TERMS 
NODES  AT  1  FT. 


II .B. — MAXIMA 


MAXIMUM 

X-COOKD 

MAXIMUM 

X-COORD 

POSITIVE 

(  IN) 

NEGATIVE 

(  IN  ) 

AXIAL  DISPLACEMENT  (IN) 

1  .  OOOE+OO 

0.00 

0. 

0.00 

LATERAL  DISPLACEMENT  (IN) 

1 .  OOOE  f 00 

0.00 

-1 .716E-02 

96.00 

ROTATION  (RAD) 

6  •  306E-04 

108.00 

-2 . 012E-02 

24.00 

AXIAL  FORCE  (P) 

0. 

0.00 

-3 . 602E+06 

0.00 

SHEAR  (P) 

2 . 077E+05 

0.00 

-4 .836E+04 

72.00 

BENDING  MOMENT  (P-IN) 

1 .709E+06 

60.00 

-5. 424E4-06 

0.00 

II .C. --REACTIONS  AT  FIXED  SUFPORTS 

X-COORD  X-REACT ION  Y-REACT ION  MOM-REACTION 

(IN)  <P>  (P>  ( P- IN  > 

0.00  3.60E+06  2.08E+05  5.42E+0& 
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PROGRAM  CBEAMC  -  ANALYSIS  OF  BEAM-COLUMNS  WITH  NONLINEAR  SUPPORTS 

hate:  06/24/80  time:  14:02:53 


III .--COMPLETE  RESULTS 

III .A. — HEADING 
EXAMPLE  3  —  CONTINUED 
AXIAL  AND  LATERAL  TERMS 
'NODES  AT  1  FT. 


III. 8. — DISPLACEMENTS  AND  INTERNAL  FORCES 


< - 1 

&ISPLACEMEN1 

S - V 

< - INTERNAL  FORCES - > 

X-COORD 

AXIAL 

LATERAL 

ROTATION 

AXIAL 

SHEAR 

MOMENT 

(  IN) 

(  IN) 

( IN) 

<RAD> 

<P> 

<p> 

( P— IN) 

0.00 

1 .000E+00 

1 . OOOE+OO 

0. 

-9. 649E+06 

2 . 077E+05 

-5 . 424E+06 

12.00 

7.350E-01 

9 . 076E-01 

-1 .402E-02 

-9. 285E+06 

1 .994E+05 

-2.966E+06 

24.00 

4.903E-01 

6.953E-01 

-2.012E-02 

-7.998E+06 

1 .655E  +  05 

-7.442E+0* 

36.00 

2 . 893E-01 

4.51  IE-01 

-1 .967E-02 

-5. 868E+06 

1 . 029E+05 

8.879E+05 

48.00 

1 .473E-01 

2 . 398E-0 1 

-1.511 E-02 

-3.575E+06 

3. 139E+04 

1 .6B9E+06 

60.00 

6.295E-02 

9.371E-02 

-9.233E-03 

-1 .763E+06 

-2.371E+04 

1 . 709E4  Oi 

72.00 

2. 167E-02 

1 .465E-02 

-4 . 205E-03 

-6. 780E+O5 

-4.836E+04 

1 .245E+06 

84.00 

S.628E-03 

-1 .476E-02 

-1 .026E-03 

-1 .913E+05 

-4.545E+04 

6.602E+05 

96.00 

9. 701E-04 

-1.716E-02 

3 . 740E-04 

-3.470E+04 

-2 . B45E+04 

2. 104E+05 

108.00 

7.709E  -05 

-1 . 036E-02 

6 . 306E-04 

-2.797E+03 

-1 . 106E+04 

-2.207E+04 

120.00 

-1 .S36L-06 

-3.864E-03 

4. 187E-04 

6.417E+01 

-4 . 502E+02 

-B.321E+04 

III.C. — FORCES  IN  DISTRIBUTED  LINEAR  SPRINGS 
DISTRIBUTED  SPRING  FORCES 


;-COORD 

AXIAL 

LATERAL 

<  IN) 

< P/ IN) 

(P/IN) 

0.00 

0. 

0. 

12.00 

-1 . 058E+03 

-1.307E+03 

24.00 

-2.824E+03 

-4.005E+03 

36.00 

-3. 749E+03 

-5.847E+03 

48.00 

-3 . 395E+03 

-5.525E+03 

60.00 

-2. 266E+03 

-3.373E+03 

72.00 

-1. 123E+03 

-7.596E+02 

84.00 

-3.971E+02 

1 .042E+03 

96.00 

-8 . 940E+0 1 

1 .582E+03 

108.00 

-8.991E+00 

1 . 208E+03 

120.00 

2.644E-01 

5.564E+02 

\  CBEAMC 

)o/24/30 


ANALYSIS  OF  BEAM-COLUMNS  WITH  NONLINEAR  SUPPORTS 

time:  14:03:43 


JHMARt  OF  RESULTS 

. --HEADING 

rL E  3  --  CONTINUED 

■IT  TERMS 


. --MAXIMA 


MAXIMUM 

X-COORD 

MAXIMUM 

X-COORD 

POSITIVE 

(IN) 

NEGATIVE 

(IN)  ; 

DISPLACEMENT  < I N  > 

0. 

0.00 

0. 

0.00 

l.  DISPLACEMENT  <IN) 

1 . 021E  +  0 1 

24.00 

-2.083E-01. 

108 . 00 

DN  .PAD' 

I . OOOE+OO 

0.00 

-2.301E-01 

48.00 

FORCE  <P> 

0. 

0.00 

0. 

0.00  1 

.  F-  » 

5.424E+06 

0.00 

-7 . 242E+05 

84.00  ; 

G  MOMENT  ( P -  I N ) 

2.380E+07 

72.00 

-2. 148E+08 

0.00 

i 


. --REAC 


IONS  AT  ^IXED  SUPPORTS 


X-COORD 

(IN) 

0.00 


X -REACT  I  ON 
(P) 

0. 


Y-RE ACT  I  ON 

<P> 

5.42E+06 


MOM-REACTION 

(P-IN) 

2 » 1 5E  +  08 


I 


l 
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PROGRAM  CBEAMC  -  ANALYSIS  OF  BEAM-COLUMNS  WITH  NONLINEAR  SUPPORTS 
BATE'.  06/24/80  TIME:  14:03:43 


III.— COMPLETE  RESULTS 

III.  A. —HEADING 
EXAMPLE  3  —  CONTINUED 
MOMENT  TERMS 


1II.B. — DISPLACEMENTS  AND  INTERNAL  FORCES 


< - 

--DISPLACEMENTS - > 

< - 

-INTERNAL  FORCES - > 

X-COORD 

AXIAL 

LATERAL 

ROTATION 

AXIAL 

SHEAR 

MOMENT 

(  IN) 

(IN) 

(IN) 

( RAD ) 

(P) 

(P) 

(P-IN) 

0.00 

0. 

0. 

1 .OOOE+OO 

0. 

5.424E+06 

-2  *  1 48E+08 

12.00 

0. 

8. 121E+00 

3.897E-01 

0. 

5.372E+06 

-1 .498E+08 

24.00 

0. 

1 .021E+01 

-5.985E-03 

0. 

4 . 950E+06 

-8.732E+07 

36.00 

0. 

8.767E+00 

-2.050E-01 

0. 

3.873E+06 

-3 . 372E+07 

48.00 

0. 

5.91 1 £+00 

-2.501E-01 

0. 

2 . 307E+06 

3 . 630E406 

60.00 

0. 

3. 136E+00 

-2.024E-01 

0. 

7 . 573E+05 

2. 173E+07 

72.00 

0. 

1 . 179E+00 

-1 .226E-01 

0. 

-3.067E+05 

2.380E+07 

84.00 

0. 

1.494E-01 

-5.290E-02 

0. 

-7.242E+05 

1 .700E+07 

96.00 

0. 

-2 . 027E-01 

-1 .062E-02 

0. 

-6. 495E+05 

8.407E+06 

108.00 

0. 

-2 . 083E-01 

6.225E-03 

0. 

-3.709E+05 

2 • 220E+06 

120.00 

0. 

-1 . 132E-01 

8.053E-03 

0. 

-1 . 192E+05 

-6.218E+05 

III.C. — FORCES  IN  DISTRIBUTED  LINEAR  SPRINGS 
DISTRIBUTED  SPRING  FORCES 
X-COORD  AXIAL  LATERAL 

(IN)  (P/IN)  (P/IN) 


0.00 

0. 

0. 

12.00 

0. 

-1 . 169E+04 

24.00 

0. 

-S . 883E+04 

36.00 

0. 

-1 . 136E+05 

48.00 

0. 

-1.362E+05 

60.00 

0. 

-1.129E+0S 

72.00 

0. 

-6. 113E+04 

84.00 

0. 

-1 .054E+04 

96.00 

0. 

1 .868E+04 

108.00 

0. 

2 . 430E+04 

120.00 

0. 

1 .630E+04 
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Deformation  Curves  for  Axial  (x)  Soil  Resistance 


*****  inFot  file  for  example  4  ***** 


E?££FLE  4  "  MULTIPLE  ANCHORED  RETAINING  WALL 
1010  EFFECTS  OF  AXIAL  FORCE  ON  BENDING  EXCLUDED 
1020  TENSION-ONLY  ANCHORS 
1030  BEAN 


1040  (PZ  38  SHEET  PILE) 

1030  0  552  29. E6  11.2  280.8 
1060  NODES  F 
1070  0  50  1 
1080  NONLINEAR  F  p 

1090  (LATERAL  SPRINGS  FOR  RIGHTSIDE  SOIL) 

tlOO  D  Y  0  2  1  1 

1  110  -0.008  0.125 

1120  -133  -1200 

1130  C  16  2  1  1 

1140  -0.008  0.125 

l  ISO  -747  -6720 

1160  C  17  2  1  1 

1170  -0.231  0.231 

1130  -303  -4303 

1190  E  46  2  1  1 

1200  -0.231  0.231 

1210  -2115  -6115 

1220  (LATERAL  SPRINGS  FOR  LEFTSIDE  SOIL) 

1230  D  Y  26  2  1  1 

1240  -0.231  0.231 

1250  2000  0 

1260  E  46  2  1  1 

1270  -0.231  0.145 

1280  3250  0 

1290  (AXIAL  SPRINGS  FOR  RIGHTSIDE  SOIL) 

1300  D  X  0  2  1  1 

1310  0  0.1 

1320  0  -61 

1330  C  16  2  1  1 

1340  0  0.1 

1350  0  -342 

1360  C  17  2  1  1 

1370  -0.1  0.1 

1380  500  -500 

1390  E  46  2  1  1 

1400  -0.1  0.1 

1410  1000  -1000 

1420  (AXIAL  SPRINGS  FOR  LEFT  SIDE  SOIL 

1430  D  X  26  2  1  1 

1440  -0.1  0.1 

1450  500  -500 

1460  E  46  2  1  1 

1470  -0.1  0.1 

1480  1000  -1000 

1490  (ANCHORS) 

1500  C  6  90  2  1  1.E3 
1510  -0.043  0 
1520  8  0 


1530  C  17  70  2  1  1.E3 
1540  -0.052  0 
1550  10  0 
1580  FINISH 


t 
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PROGRAM  CBEAMC  -  ANALYSIS  OP  BEAM-COLUMNS  UITH  NONLINEAR  SUPPORTS 
DATE:  06/24/80  TIME:  19:36:42 

I. --INPUT  DATA 

1. --HEAPING 

EXAMPLE  4  —  MULTIPLE  ANCHORED  RETAINING  UALL 
EFFECTS  OF  AXIAL  FORCE  ON  BENDING  EXCLUDED 
TENSION-ONLY  ANCHORS 


2. --BEAM  CROSS  SECTION  DATA 


X-COORDINATE 

MODULUS  OF 

< - ST 

---SE-TION 
ART - 

PROPERTIES 

STOP - 

START 

STOP 

ELASTICITY 

AREA 

INERT. A 

AREA 

INERTIA 

(  IN) 

(  IN) 

(PSD 

(SOIN) 

< IN**4) 

( SQIN ) 

( IN$*4  ) 

0.00 

552.00 

2.90E+07 

it  *20 

280.80 

11.20 

O 

CD 

O 

CO 

r< 

3. --NODE  SPACING  DATA 

X-COORDINATE  MAXIMUM  NODE 

START  STOP  SPACING 

(FT)  (FT)  (FT ) 

0.00  50.00  1.00 


4. --LOAD  DATA 
NONE 


5. --FIXED  SUPPORT  DATA 
NONE 


6. --LINEAR  SPRING  DATA 
NONE 


--NONLINEAR  SPRING  DATA 


7. A. --NONLINEAR  CONCENTRATED  SPRINGS 


SPRING 

CURVE  COORD. 

MULTIPLIERS 

X-COORD 

ANGLE 

DISPLACEMENT 

FORCE 

(FT) 

(DEG) 

<FT  ) 

<P  ) 

6.00 

90.00 

1 .00E-03 

1.00E+03 

DISP.  COORDS 

.  -43.00  0.00 

FORCE  COORDS 

.  8.00  0.00 

SPRING 

CURVE  COORD. 

MULTIPLIERS 

X-COORD 

ANGLE 

DISPLACEMENT 

FORCE 

<  FT  > 

(DEG) 

(FT  > 

<P) 

17.00 

70.00 

1 .00E-03 

1 .00E  +  03 

DISP.  COORDS 

.  -52.00  0.00 

FORCE  COORDS 

♦  10.00  0.00 

7.P.--N0NL 

INEAR  DISTRIBUTED 

SPRINGS 

SPRING 

SPRING 

CURVE  COORD. 

MULTIPLIERS 

DISTRIBUTION  DIRECTION 

X-COORD 

DISPLACEMENT 

FORCE 

(FT) 

(FT  ) 

(P) 

STARTS  Y 

0.00 

1 .OOE-02 

1 . 00E+02 

DISP. 

COORDS.  -.80  12.30 

FORCE 

COORDS.  -1.33  -12.00 

SPRING 

SPRING 

CURVE  COORD. 

MULTIPLIERS 

DISTRIBUTION  DIRECTION 

X-COORD 

DISPLACEMENT 

FORCE 

(FT) 

(FT) 

(P> 

CONTINUES  Y 

16.00 

1 .00E-02 

1 . 00E+02 

DISP. 

COORDS.  -.80  12.50 

FORCE 

COORDS.  -7.47  -67.20 

SPRING 

SPRING 

CURVE  COORD. 

MULTIPLIERS 

DISTRIBUTION  DIRECTION 

X-COORD 

DISPLACEMENT 

FORCE 

(FT) 

(FT) 

(P) 

CONTINUES  Y 

17.00 

1 .00E-O2 

1  .00E  +  02 

DISP. 

COORDS.  -23.10  23.10 

FORCE 

COORDS.  -3.03  -43.03 

SPRING 

SPRING 

CURVE  COORD. 

MULTIPLIERS 

DISTRIBUTION  DIRECTION 

X-COORD 

DISPLACEMENT 

FORCE 

(FT) 

(FT) 

(P> 

ENDS  Y 

46.00 

1  .00E-02 

1 .00E+02 

DISP.  COORDS.  -23 . 10  23.10 
FORCE  COORDS.  -21.13  -61.13 


SPRING 

SPRING 

CURVE  COORD. 

MULTIPLIERS 

DISTRIBUTION  DIRECTION 

X-COORD 

DISPLACEMENT 

FORCE 

(FT  > 

(FT  ) 

(P) 

STARTS  Y 

26.00 

1  .OOE-02 

1 . 00E  +  02 

DISP.  COORDS.  -23.10  23.10 

FORCE  COORDS.  20.00  0.00 

SPRING 

SPRING 

CURVE  COORD. 

MULTIPLIERS 

DISTRIBUTION  DIRECTION 

X-COORD 

DISPLACEMENT 

FORCE 

(FT  ) 

(  FT  ) 

(P) 

ENDS  Y 

46.00 

1 .00E-02 

1 .00E  +  02 

DISP.  COORDS.  -23.10  14.30 

FORCE  COORDS.  32.50  0.00 

SPRING 

SPRING 

CURVE  COORD. 

MULTIPLIERS 

DISTRIBUTION  DIRECTION 

X-COORD 

DISPLACEMENT 

FORCE 

<  FT  ) 

(FT) 

(P  > 

STARTS  X 

0.00 

1  .O0E-03 

1  .OOE+OO 

DISP.  COORDS.  0.00  100.00 

FORCE  COORDS.  0.00  -61.00 

SPRING 

SPRING 

CURVE  COORD. 

MULTIPLIERS 

DISTRIBUTION  DIRECTION 

X-COORD 

DISPLACEMENT 

FORCE 

(FT  ) 

(FT) 

(P) 

CONTINUES  X 

16.00 

1 .00E-03 

1 .00E+01 

DISP.  COORDS.  0.00  100.00 

FORCE  COORDS.  0.00  -34.20 

SPRING 

SPRING 

CURVE  COORD. 

MULTIPLIERS 

DISTRIBUTION  DIRECTION 

X-COORD 

DISPLACEMENT 

FORCE 

(FT) 

(FT) 

(P) 

CONTINUES  X 

17.00 

1 .00E-O3 

1  .OOE  +  Ol 

DISP.  COORDS.  -100.00  100.00 

FORCE  COORDS.  50.00  -50.00 

SPRING 

SPRING 

CURVE  COORD. 

MULTIPLIERS 

DISTRIBUTION  DIRECTION 

X-COORD 

DISPLACEMENT 

FORCE 

(FT) 

(FT) 

(P) 

ENDS  X 

46.00 

1  .  OOE-03 

1 .00E  +  02 

DISP.  COORDS.  -100.00  100.00 

FORCE  COORDS.  10.00  -10.00 

SPRING 

SPRING 

CURVE  COORD. 

MULTIPLIERS 

DISTRIBUTION  DIRECTION 

X-COORD 

DISPLACEMENT 

FORCE 

(FT) 

(FT) 

(P) 

STARTS  X 

26.00 

1 .00E-03 

1 . OOE+Ol 

DISP.  COORDS.  -100.00  100.00 

FORCE  COORDS.  50.00  -50.00 

SPRING 

SPRING 

CURVE  COORD. 

MULTIPLIERS 

DISTRIBUTION  DIRECTION 

X-COORD 

DISPLACEMENT 

FORCE 

(FT) 

(FT) 

(P) 

ENDS  X 

46.00 

1. 00E-03 

1 .00E+02 

DISP.  COORDS.  -100.00  100.00 
FORCE  COORDS.  10.00  -10.00 
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PROGRAM  CBEAMC 
HATE!  06/24/80 


ANALYSIS  OF  BEAM-COLUMNS  WITH  NONLINEAR  SUPPORTS 

time:  i ? : 3 7 : 2 i 


1 1. --SUMMARY  OF  RESULTS 


II  .A. --HEADING 

EXAMPLE  4  --  MULTIPLE  ANCHORED  RETAINING  UALL 
EFFECTS  OF  AXIAL  FORCE  ON  BENDING  EXCLUDED 
TENSION-ONLY  ANCHORS 


I  I  ,B. --MAXIMA 


MAXIMUM 


X -COORD 


AXIAL  DISPLACEMENT  (FT) 
LATERAL  DISPLACEMENT  (FT) 
ROTATION  (RAD) 

AXIAL  FORCE  <P> 

SHEAR  (P) 

BENDING  MOMENT  (P-FT) 


POSITIVE 
8.584E-03 
1.775E-02 
5.487E-05 
3. 128E  +  02 
7 « 412E+03 
3 . 6 1 8E+ 04 


II  .C  .--FORCES  IN  CONCENTRATED  NONLINEAR  SPRINGS 
X-COORD  ANGLE  DEFORMATION 


(FT) 

6.00 

17.00 


(DEG) 

90.00 

70.00 


(FT  ) 

-3 . 452E-02 
-1  . 145E-01 


MAXIMUM 


X-COORD 


(FT  ) 
17.00 
0.00 
46.00 
17.00 
17.00 
23.00 


NEGATIVE 


.723E-01 
.  437E-04 
.  107E  +  03 
.  155E  +  03 
.  6  38E  +  0  3 


(FT  ) 
0.00 
32.00 
7.00 
17.00 
26.00 
6.00 


FORCE 

(P) 

6.423E+03 
1 • 000E+04 


PROGRAM  CBEAMC  -  ANALYSIS  OF  BEAM-COLUMNS  WITH  NONLINEAR  SUPPORTS  BM 
DATE!  06/24/80  TIME!19!37!21 

III .--COMPLETE  RESULTS 

III .A. --HEADING 

EXAMPLE  4  --  MULTIPLE  ANCHORED  RETAINING  WALL 
EFFECTS  OF  AXIAL  FORCE  ON  BENDING  EXCLUDED 
TENSION-ONLY  ANCHORS 


III. B. --DISPLACEMENTS  AND  INTERNAL  FORCES 
- UISPl  ACEMENTS - > 


- INTERNAL  FORCES- 


COORD 

AXIAl 

LATERAL 

ROTATION 

AXIAL 

SHEAR 

MOMENT 

(FT) 

(FT) 

(FT) 

(  RAD) 

(P  ) 

(P  ) 

(F-FT  ) 

0.00 

8.578E-03 

1 .775E-02 

-8 . 660E-03 

0. 

0. 

0. 

1  .00 

8 . 57 BE -0  3 

9.087E-03 

-8 .661E-03 

5.986E+00 

-3 • 470E+02 

-I  . 728E  +  02 

2.00 

8 .578003 

4 . 236E  -04 

-8.668003 

1 .348E  +  01 

-6 . 824E+02 

-6 . 902E+02 

3.00 

8.578003 

-8.253E-03 

-8 . 687E-03 

2.248E+01 

-9.463E+02 

-1 .510E  +  03 

4.00 

8.578003 

-1 .696  0  02 

-8 . 723E-03 

3.29BE+01 

-1 .214E+03 

-2 . 587E+03 

5.00 

8 . 578003 

-2.571E-0? 

- 8 . 780E-03 

4 . 499E+01 

-1 .519E  +  03 

-3.950E+03 

6.00 

8 .578003 

3.452E-02 

•8 .8o4E-03 

5.851E+01 

-1 .863E+03 

-5 .638E+03 

6.00 

8.5781 -03 

-3.452E  02 

-0.864E -03 

5. 05 IE +01 

4 .560E  +  03 

-5.638E+03 

7.00 

8.579E -03 

-4 . 343E-02 

-8.925E-03 

7. 354E+01 

4. 177E  +  03 

-1 . 266E+03 

8.00 

8 . 579003 

-5.235E -02 

8.912E-03 

9.007E+01 

3.757E+03 

2 . 704E+03 

9,00 

8.579  0  03 

6. 123F.-02 

-8. 832 E -03 

1  » 081 E  +  0? 

3  *  297E+03 

6 . 234E+03 

10.00 

8.579E *03 

-6.999E- 02 

- 8 . 694C-03 

1 .277L+02 

2 .800E  +  03 

9.286E+03 

11.00 

8 . 580003 

-7.860E  02 

-8 .506F-03 

1 . 487E+02 

2 . 264E+03 

1  .  182E+04 

12 . 00 

8.580E-03 

-8. 6991 -02 

-8 . 279L-03 

1 .713E  +  0? 

1 .690E  +  03 

1 .380E  +  04 

13.00 

8 .531F  03 

-9.5151  -02 

-8 . 022E-03 

1 .953F+0? 

1 .077E+03 

1 .519E  +  04 

14.00 

8  •  582003 

-1 .030E-01 

-7.746E-03 

2  » 209E+02 

4.259E+02 

1 .594E+04 

15.00 

3.582E--0J 

-1  . 106E-01 

-7. 462E-03 

2.480E+02 

-2.63SE+02 

l  .603E  +  04 

16.00 

8  <j58 3E  -  03 

-1  . 180E  -01 

-7. 183E-03 

2. 766E+02 

-9.91 3E  +  02 

1 . 540E+04 

17.00 

8.584003 

-1 .250E-01 

-6.922E-03 

3. 128E  +  02 

-1 .985E+03 

1 .395E  +  04 

17,00 

8.594E-03 

-1 . 250001 

-6 . 922E-03 

-3. 107E+03 

7.412E+03 

1 .395E+04 

18.00 

8.575E-03 

-1 .318E-01 

-6.613E  -03 

-3. 064E+03 

6. 189E+03 

2.075E+04 

19.00 

8.565E-03 

-1 .382E-01 

-6. 195E-03 

-3.019E+03 

4 .962E  +  03 

2.633E+04 

20.00 

8.556E  03 

-1 . 441E  -01 

-5.689E-03 

-2.972E+03 

3.725E+03 

3  •  067E+04 

21 .00 

8.547E-03 

-1 .496001 

-5. 110E-03 

-2 . 924E+03 

2.475E+03 

3 . 378E+C4 

22.00 

8.538E-03 

-1 .544001 

-4.502E  -03 

-2.875E+03 

1 .207E  +  03 

3 . 562E+04 

23.00 

8  •  5  2  9  E  -  0  3 

-1 . 5860  01 

-3.865E-03 

-2.824E+03 

-8  »  444E  +  01 

3.618E+04 

24.00 

8. 52 IE -03 

-1 .621001 

-3.230F-03 

-2 . 772E+03 

-1 .405E+03 

3.544E+04 

25.00 

8.512E- 03 

-1 .650001 

-2 . 620E-03 

-2. 71BE+03 

-2. 760E+03 

3.336E+04 

26.00 

8.50  40  03 

-1,6741  -01 

-2.058E-03 

-2.663E+03 

-4 . 155E+03 

2.991E+04 

27.00 

8.496E-03 

-l .692001 

-1 .565E-03 

-2.563E+03 

-3 . 84 1 E  +  03 

2.591E+04 

28.00 

8.488003 

-1  .  705E-01 

-1  .  140E-03 

-2.460E+03 

-  3 . 520E+03 

2. 223E  +  04 

2*. 00 

8.4R1O02 

-1 . 715E  01 

-7 . 7 7 3 E - 0 4 

353E+03 

-3. 197E  +  03 

1 .887E+04 

30.00 

8.474F-03 

-1 . 72  IF -01 

-4 . 7 1  OF -04 

-2 . 243F+03 

-2  »  876E  +  03 

1 .583E+04 

3 1  .  Oo 

8 . 46?[.  -03 

-1 . 7241-01 

-2. 155E-04 

-2. 129E+03 

-2 . 561E+03 

1 .  3 1 1 E  +  04 

32.00 

3. 460F  -03 

-1 . 72500  1 

-5.342006 

-2.01 1E  +  03 

-2, 255E+03 

1  »  0  7 1 E  +  0  4 

33.00 

3  *  454003 

-.1 .7740  01 

1 .6491.-04 

-1 .890E+03 

-1 .  962E  +  0 3 

8 . 599E+03 

34.00 

8 . 449F-03 

-1 .722001 

3.005E-04 

-1 . 766E+03 

-1 .684E+03 

6 . 77  7E  +0  3 

35.00 

8. 444E  03 

-1  .  719001 

4 . 062E-04 

-1 • 63BE+03 

-1 . 423E+03 

5  •  225E+03 

36.00 

8 . 4  3VE -03 

-1  . 714E  -01 

4.867E--04 

-1 . 507E+03 

-1  .  181E  +  03 

3.924E+03 

37.00 

.  *  34003 

-  1 . 709E-01 

5. 464E-04 

-1 .372E  +  03 

-9.593E+02 

2 . 855E+03 

38.00 

8.  30E -03 

-1  .  703E-01 

5 .890E-04 

-1 .233E+03 

-7.587E+02 

1 .998E+03 

39.00 

8. 427! -03 

-1 . 6  9  7  E - 0 1 

6. 1820  0  4 

-1 .091E  +  03 

-5.803E+02 

1 . 33 1 E  +03 

40 . 00 

8.424E-03 

-1.69  IE -01 

6. 370E--04 

-9.461Ef02 

-4 . 248E+02 

8  «  30 l E  +  02 

41.00 

8 . 421E-03 

-  1 . 6H4C-01 

6,484004 

-7 »  9 72C  +0 2 

-2.928C+02 

4. 733E+02 

42.00 

3. 4( 9003 

1 . 679E -01 

6.545E-04 

-6.448E+02 

-1 .849E+02 

2.364E+02 

43.00 

0 . 41 7T-03 

-1 .671001 

6 . 57  3004 

-4.889E+02 

-1 .013E+02 

9.538E+01 

44.00 

8 . 416E -03 

-1 .665E-01 

6 .583F- 04 

-3.295E+02 

-4.244E+01 

2.557E+01 

45.00 

8 . 4 1 50  03 

-1 . 65CE-01 

6.585E-04 

-1 . 665E+02 

-8 .572E+00 

2. 166E+00 

46.00 

8 .415E-03 

-1 .652E-01 

6.585E-04 

0. 

0. 

0. 

Ill .C. --FORCES  IN  DISTRIBUTED  NONLINEAR  SPRINGS 
DISTRIBUTED  SPRING  FORCES 


COORD 

AXIAL 

LATERAL 

(FT; 

(P/FT ) 

(P/FT) 

0.00 

-  5 . 2  3  3  E  +  0  0 

-3. 396E  +  02 

1 .00 

-  6 . 739E+00 

-3.4“?9E  +  02 

2.00 

-8 . 246E+00 

-3. 162E+02 

3.00 

-9. 75 2E +00 

-2. 481 E +  02 

4 . 00 

1 . 126E+01 

-2 . 865E+02 

5.00 

■1 .277E+01 

-3.249E+02 

6. 00 

-1 .427E+01 

-3 . 633E+02 

7 . 00 

-1 .578E+01 

-4 . 016E+02 

8  .  00 

~1 • 729E+01 

-4 . 400E+02 

9.00 

-1 .879E+01 

-4 . 784E+02 

1.0.00 

~2 . 030E+01 

-5. 168E+02 

1  1 . 00 

--2. 181  F +  01 

-5.551E+02 

12.00 

-2.3321+01 

-5 . 935E  +  02 

13.00 

-2.403E+01 

-6.319E+02 

14.00 

-2.633E+01 

-6.703E+02 

15.00 

-2.784E+01 

-7.086E  +02 

16.00 

-  2 . 9  3  5  E  +  0 1 

-  7 . 470E  +  02 

17.00 

4.292E+01 

-1 .221E+03 

18.00 

-4. 43 5 E +  01 

-1 .224E+03 

19.00 

-4 . 578E  +  01 

-1 . 231E  +  03 

20 . 00 

-4.721E+01 

-1 . 242E+03 

21 .00 

-4 . 863E+01 

-1 . 258E+03 

’2.00 

-  5 . 0  0  5  E  +  0 1 

-  1 .279E  +  03 

23.00 

-5.147E+01 

-1 .305E+03 

24.00 

-  5 . 2  8  9  F  +  0  1 

-1 .337E+03 

25.00 

-5.430E+01 

-  1 . 3/4E  +  03 

7  A  .  00 

-5 . 57 IE +  0.1 

-1 .416E+03 

26.00 

-9.823E+01 

3.081E+02 

27.00 

-1 . 01 7E  +  02 

3 . 1  82F  +  02 

28.00 

-1 .052E  +  02 

3 . 230E+02 

29.00 

-1 . 087E+02 

3 , 229E+02 

30.00 

-1 . 122E  +  02 

3 . 187E+02 

31 . 00 

1  . 157E  +  02 

3. 108E+02 

32.00 

-1  .  192E  +  02 

2.997E+02 

33.00 

-  1 . 2  2  7  E  +  0  2 

2.859E+02 

34.00 

-1 .261E  +  02 

2. 698E+02 

35.00 

-1 . 2  9  6  E  +  0  ? 

2. 518E+02 

36.00 

-1 . 331E  +  02 

2 . 3  2  3  E  +  0  2 

37.00 

-1 . 366E+02 

2.115E+02 

38.00 

-1 . 40  IF. +  02 

1 . 897E+02 

39.00 

-  1  . 436E  +  02 

l .67] E+02 

40 . 00 

-1  .  4  7  1 E  +  0  2 

1 . 438E+02 

41.00 

-1 . 5  0  6  E  +  0 2 

1 . 200E+02 

42 . 00 

-1.  .54  If  902 

9 . 582E+01 

43.00 

-1.  ,  5  7  7  £  -1  0  2 

7 . 1 28E+01 

44.0c 

-1 .6121+02 

4 . 642E+01 

45.00 

-1 . 6  4  7  E 1 0 2 

2. 127E+01 

'*6.00 

-1  •  603E  +  02 

-4 . 171E+00 
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In  accordance  with  letter  from  DAEN-RDC,  DAEN-ASI  dated 
22  July  1977,  Subject:  Facsimile  Catalog  Cards  for 
Laboratory  Technical  Publications,  a  facsimile  catalog 
card  in  Library  of  Congress  MARC  format  is  reproduced 
below. 


Dawkins,  William  P. 

User's  guide:  Computer  program  for  analysis  of 
beam-column  structures  with  nonlinear  supports  (CBEAMC)  /  by 
William  P.  Dawkins.  --  Vicksburg,  Miss.  :  U.S.  Army 
Engineer  Waterways  Experiment  Station  ;  Springfield,  Va.  : 
available  from  NTIS,  1982. 

90  p.  in  various  pagings  ;  ill.  ;  27  cm.  --  (Instruction 
report  ;  K-82-6) 

Cover  title. 

"June  1982." 

Final  report. 

"Prepared  for  Office,  Chief  of  Engineers,  U.S.  Army." 
"Monitored  by  Automatic  Data  Processing  Center,  U.S. 

Army  Engineer  Waterways  Experiment  Station." 

"A  report  under  the  Computer-Aided  Structural 
Engineering  (CASE)  Project." 

Bibliography:  p.  36. 

1.  CBEAMC  (Computer  program).  2.  Computer  programs. 

3.  Mathematical  models.  4.  Structural  frames--Models . 


Dawkins,  William  P. 

User's  guide:  Computer  program  for  analysis  :  ...  1982. 

(Card  2) 

5.  Structures,  Theory  of.  I.  United  States.  Army.  Corps 
of  Engineers.  Office  of  the  Chief  of  Engineers.  II.  U.S. 
Army  Engineer  Waterways  Experiment  Station.  Automatic  Data 
Processing  Center.  III.  Title  IV.  Series:  Instruction 
report  (U.S.  Army  Engineer  Waterways  Experiment  Station)  ; 
K-82-6. 
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